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ABSTRACT 


Model-generated data and FGGE analyses are used to 
calculate quasi-Lagrangian budgets of mass, vorticity, heat 
and moisture following disturbances that form within podas 
air streams. These polar lows grow primarily through basic 
baroclinic instability processes and exhibit many features 
of larger maritime extratropical cyclones. Polar lows that 
originate on the equatorward side of a nearly straight 
upper-level jet are contrasted with lows that form on the 
poleward side of the jet and have considerable  mid- 
tropospheric positive vorticity advection at formation time. 
The absence of favorable coupling to a jet stream was the 
missing factor in two model-generated polar lows Enat faired 
to develop. Although the vorticity balance 1s initial 
different for the two types of polar lows, the vortice 
budgets during later stages are similar. The heat budget and 
the thickness tendency equation demonstrate that the self- 
development process that is present in larger maritime 


cyclones is also important for polar low intensif o sr n un 
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I. INTRODUCTION 


599 POLAR AIR STREAM CYCLOGENESIS 


The existence of cyclogenesis in polar air streams 
upstream or poleward of the main polar frontal cloud band 
has long been recognized by operational meteorologists. 
Such disturbances form primarily over the ocean in winter 
(Reed, 1979), and are of relatively small size (500-1500 km) 


compared with the typical extratropical cyclone. Over the 
ocean, their inciplent stage is sometimes characterized by 
areas of locally enhanced cumulonimbus activity. Mis 


convective activity may later become organized into a comma- 
кәс cioud pattern (Reed, 1979; Mullen, 1983), or acquire 
£- Pres similar to tropical storms (Rasmussen, 1979), The 
disturbance is often accompanied by appreciable convective 
EEPcipitation, and occasionally by severe veather 


1— on 1968: Suttie, 1970; Lyall, 1972; Seaman et al., 


1982). Cyclogenesis is Sometimes vigorous with considerable 
areal expansion of the disturbance (e.g., Cook, |l. 
Mullen, 1983). The terms polar low, polar air depression, 
pur trough, cold trough Şoaırmass trough, comma cloud, 
Epor disturbance, cold vortex, Arctic instability low, and 
Arctic Bomb have been applied to this phenomenon. The term 


B— “lov will here refer to any cyclonic disturbance 
Spreinating within a polar air stream. 

Polar low occurrence is more widespread and frequent 
than the relatively sparse treatment in the literature would 
muy. Weber and Wilderotter (1951) state that Polar- air 
stream cyclogenesis is the most frequent type observed over 
acne North Atlantic and North Pacific Oceans, and the 


12550 American continent. 1205505107” (017770) shovs that 


иә 


between 20% and 50% of the annual San Francisco precipita- 
tion occurs in association with non-frontal air mass distur- 
bances akin to the polar low. Cyclogenesis away from pre- 
existing cloud bands accounts for about 504 of alm Sot niema 
Hemisphere winter cyclogenesis (Streten and Troup, 1973). 
Clark (1977) estimates that polar lows are involved in two- 
thirds of all winter cyclogenesis events in the Southern 
Hemisphere, where development is frequently of the so-called 
"instant occlusion mode (Carleton hm... An instan 
occlusion forms when a polar low (evidenced by organized 
convection, usually in a comma formation) advances on an 
incipient frontal wave development. The resulting cloud 


vortex rapidly develops the signature of a mature occluded 


system (Anderson et al., 1969). Sometimes the comma cloud 
and frontal wave merge. In other cases, the comma cloud 
develops into the occlusion, distinct from the main frome 


(Locatelli et wi. oa) 

It has been suggested that two different types of poem 
lows exist (Rasmussen, 1979, 1981; Locatelli et 4123) 71702 
Sarde. and Warner, 1933). Satellite imagery indicates thas 
polar lows exhibit a variety of cloud signatures (Forbes and 
Lottes, 1982). Some bear the hallmarks of organized convec- 
tive tropical cloud systems. Okland (IPP) and Rasmuss m 
(1979) applied thermal instability theory to these systent 
and proposed that heating and cumulus convection play essen- 
tial roles in their ~deve lopment These convective distur- 
Dances form within the polar air mass well to Che nortiu r. 
the polar jet. 

The second type of polar low 1S more baroclın 22772 
nature. In a case study ofa polar low that cross Ди. 
British Isles, Harrold and Browning (1969) more ии 
enhanced baroclinity was evident at early stages. In addi- 
tion, the area of significant precipitation vas pProduce 22 


Stable large-scale ascent rather than by small-scale 


LA 


Convection. Suluğubaroclinic disturbances are found to 
develop in the vicinity of the polar jet, and to exhibit 
cloud features typical of larger extratropical cyclones 
— :7177:7: Mullen, 1979, 1983; Locatelli et al., 1982, 
Cook, 1983). Because these baroclinic disturbances are 
common in the Pacific Ocean, they are frequently called 


Eocific polar lows. 


1205 PURPOSE OF STUDY 


The purpose of this study is to investigate the reasons 
Roche formation and growth of the Pacific polar low. rt 


has been proposed that baroclinic instability modified by 


diabatic processes, is responsible for the growth of the 
Ы Ce polar low (Sardie and Warner, 1983; Hodur, 1984). 
According to Reed (1979), there is little evidence to 


support the view that the Pacific polar low is a purely 
convective, barotropic phenomenon. However, the larger-scale 
Smeal zation of convection which is frequently observed in 
conjunction with polar lows may contribute to development. 
TRUS, T U IS inportant to Consider the effects of moist 
processes on development. 

The early studies of baroclinic instability (Charney, 
1947; Eady, 1949) showed that the longer wavelengths (>4000 
km) are the most unstable, and observations verify that this 
scale is most prevalent in mid-latitudes. The small scales 
Ae olar lows compared with typical mid-latitude baroclinic 
disturbances may be explained by the small static stabili- 
ties which exist at low levels. This explanation was tested 
£7 ssfield (1974) and Duncan (1977) who applied model 
Eeeremlations to observational data from polar low cases, and 
found predicted wavelengths close to those observed. 

25 her Support for baroclinic instability being the 


primary mechanism for polar low development is offered in 


185 


recent theoretical studies. Gall (1976a) and Staley and 
Gall (197/) used linear models to shov that flovs having a 
Richardson number significantly lover near the surface than 
aloft favor higher grovth rates for short vaves. Since the 
Richardson number measures the ratio of static stability to 
the square of the vertical vind shear (i.e., baroclinity), 
this theory suggests that the small horizontal scales of 
polar lovs result from either Treduced static stability or 
enhanced baroclinity (or both) near the surface as compared 
E daloTt: 

Linear baroclinic cheoryrusama uniform Zonal rm 
predicts that short baroclinic waves are confined to near 
the surface (Staley and Gall, 1977; Mansfield; 1974) a 
is at variance with polar low observations (Harrold 2 
Browning, 1969; Reed, 1979; Mullen, 19/79) Which Te Ses 
Significant geopotential perturbations near the tropopause. 
Tokioka (1973), Gall (1976b) and Hodur (1984) found €h::22 
inclusion of latent heat resulted in a deeper disturbance, a 
shift of the maximum growth rates to shorter waves, and a 
general increase in growth rate. small-scale cyclossmgw 
development occasionally occurs over the winter continent misi 
the absence of significant moist convection or surta 
Eluses: (Mullen AS I Additional theory is needed to 
explain the small scales and vertical depth of these drier 
drocıurlancees: 

Reed (1979) suggested that the comparison of pola R 
structure with results from linear baroclinic theory ise 
not be valid because polar lows frequently develop on 
already perturbed basic states. Frederiksen (1979) showed 
that the presence of basic state planetary waves causes the 
maximum amplitude for linear baroclinic waves to "voce mm 
the upper troposphere as opposed to purely zonal large-scale 
States which generally exhibit maximum amplitudes near the 


str bade” Blumen (1980) showed that nonlinear interact 


16 


between growing long and short baroclinic vaves can produce 
a deep disturbance of short horizontal scale. 

A goal of this research is to clarify the relative 
importance of upper and lov-level forcing in polar lov 
formation. The studies discussed earlier indicate that the 
occurrence of an upper-level perturbation is a consequence 
of baroclinic development and is related to the effects of 
latent heating and nonlinear processes. The presence of an 
upper-level disturbance at later stages of polar lov devel- 
opment is apparent from the associated comma-shaped cloud 
pattern. Such features are frequently referred to as upper- 
level vorticity maxima by forecasters (Reed, 1979). However, 
the role of upper-level features in early stages of polar 
low formation is not clear. The polar lows described by Reed 
(1979) and Mullen (1979) form on the poleward side of the 
upper-level wind maximum under cyclonic upper-level flow. In 
Semenast, the polar low of Cook (1983) initially formed 
beneath anticyclonic flow on the equatorward side of the 
HET yet. It is apparent that in some circumstances, polar 
lows form as a result of low-level baroclinic processes 
5: msfield, 1974; Harold oran 1769) whi lemin 
others, upper-level forcing is present (Reed, 1979). 

we aditional goal of this work is to clarify the role 
277 "ace fluxes in polar low formation and development. 
Polar lows form mostly over the oceans in winter (Reed, 
EN when surface heating is largest. Mullen (1979) showed 
that large air-sea temperature differences accompany Pacific 
polar lows, especially at genesis. Ti od i on, numerical 
Simulations have shown that latent heat release and surface 
fluxes are essential for polar low development (Sandgathe, 
meee, and for realistic structure and growth rates (Hodur, 
BI On the other hand, Mansfield (1974) found that 
surface fluxes hinder polar low development by adding heat 


in antiphase to the growing temperature wave. Danard and 
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Ellenton (1980) suggested that heating in the cold sector of 
a developing cyclone destroys available potential energy. 
These authors and Bosart (1981) suggested that the role of 
the surface fluxes is to establish a favorable low-level 


environment for polar low development. 


CS Y PONHES ES 


The main hypothesis of this study is that Pacific-type 
polar lows form as a result of moist baroclinic instab11177 
As such, they exhibit many typical features of extratropiaam 
cyclones. It is proposed that the baroclinic polar low =i m 
form initially under straight upper-level flow, or alte 
nately, downstream from a significant pre-existing upper- 
level trough. Subsequent development may be described in 
terms of the  "self-development'" mechanism. The upward 
sensible heat flux maintains the reduced low-level static 
Stability which is required for continued development of the 
disturbance. The surface moisture flux ensures an abundant 
supply of moisture for condensation which consideraba 
enhances cyclogenesis by the enhancement of vertical ascent 


and generation of available potential energy. 


To test this hypothesis, a diagnostic analys P 
several model-generated and observed polar Hows 1s 
performed. The simulated polar lows formed realistically aims 
the polar air stream of a mature cyclone. Therefore, the 


effects of nonlinear interactions and a perturbed Лав 
scale state, previously shown to be important, are presenta 
from genesis. By contrast, many modelling studies Start uM 
a uniform zonal initial state. hHodel outpullrorms 4544605855 
data set accurate to computer roundoff, with much improve 
Spatial and temporal resolution over observation BE G> m 
Adequate resolution is crucial for such a Study on soca 
lows. Conclusions drawn from study of the mode reac oe ME 


be validated by similar analysis of observed poli mEnE 
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5007000077 formation and evolution is interpreted ina 
vorticity and thermodynamic balance context. 
Storm-following budgets of vorticity and heat are used to 
diagnose important physical processes and are related via 
hydrostatic arguments. A moisture budget aids interpretation 
of diabatic terms in the heat budget. The analysis is 
complete, because the vorticity and thermodynamic equations 
emt ain no assumptions or approximations. 

The study consists of four sections. Chapter II 
describes the quasi-Lagrangian diagnostic analysis applied 
to examples of polar low development as simulated by the 
United States Navy Operational Global Atmospheric Prediction 
System (NOGAPS) model (Sandgathe, 1981), and occurring in 
the First GARP Global Experiment (FGGE) data set. Chapter 
III describes examples of polar lows from both the FGGE data 
and the NOGAPS simulations which develop under straight 
upper-level flow. Polar lows which. form in the presence of a 
pre-existing upper-level trough will be analyzed in Chapter 
IV. Finally, these developing disturbances will be compared 
with non-developing disturbances to understand why some 
disturbances grow and others do not (Chapter V). Numerical 
details of the quasi-Lagrangian diagnostics are presented in 


the Appendix. 


j) 


II. NUMERICAPTPRGGSEDURE 


This study of polar lows is based on both the results 
numerical simulations and observational FGGE data. The 
rationale for using model-generated data is that the distri- 
bution of conventional observations over the oceans is inad- 
equate to resolve polar lov disturbances. Model output 
forms a complete data set; that is, fields of meteorological 
variables which are consistent solutions to the primitive 
equations are available at every gridpoint to computer roun- 
dorf accuracy: The atmospheric model (NOGAPS) includes a 
Sophisticated parameterization of boundary layer and 
diabatic processes, vhich have been shovn to be important 
for polar lov development. Such a complete and accurate 
specification of physical. processes cannot be achieved in an 
observational study. 

In observational studies, many important processes ame 
diagnosed indirectly. Direct measurements of latent heat 
release, surface fluxes, radiative effects and friction 25. 
seldom available. Rather, these processes are usually esti- 
mated from the observed fields of wind, temperature and 
moisture. For example, diabatic effects are usually calcu- 
lated as a heat budget residual using an estimate of 
vertical velocity derived from the horizontal vind fields. 
Surface fluxes of heat, moisture and momentum (friction) are 
diagnosed via empirical relationships. In the NOGAPS model, 
these internal processes are calculated directly ino a 
arly consistent, fashion. and are available at every grid- 
point to computer roundorl acenına.. Thus, it is possibim 
to study quantitatively the effect of these process p 
model storm development (e.g. Sandgathe, 1981,  Winninghoftf 
and Elsberry, 1983; Bosse, 1984: Liou and ElsSber 75 72 
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Another important advantage to using model-generated 
data is the increased time resolution. Polar lovs are typi- 
cally short-lived so the coarse (12 h) time resolution of 
most observational analyses is often inadequate to fully 
describe rapid changes. The increased precision of time 
derivatives coupled with the dynamic consistency of the 
model-generated fields should yield reliable diagnoses of 
model storm processes. The residuals of budget-type calcula- 
tions in which a time derivative is related to the sum of 


several processes should provide an estimate of the integ- 


1215077 of the diagnostic routines. Error also arises from 
numerical approximations within the model (e.g. finite 
differences), so the residual only provides an upper bound 


on numerical uncertainty associated with the diagnostic 
foutines. 

It should be emphasized that the fidelity of the model 
Simulations in describing natural processes is strongly 
model dependent. The model is at best an incomplete repre- 
sentation of the atmosphere, and results are affected by 


numerical solution techniques. Physical processes are param- 


eterized in the model. Although computed to computer round- 
accuracy, they are only approximations Eo. troe 
atmospheric physical processes. Therefore, the study “ot 


model output will only serve as a guide to the study of 
observational data. Conclusions based on the analysis of the 
model are valıdated only vhen corroborated by similar 
results from the observational data. FGGE analyses provide 
the best oppor tunity for such a study of maritime 
cyclogenesis. 

ine tollowing sections describe the NOGAPS model and 
Bret details of the simulations, the FGGE data set and the 


basis of the quasi-Lagrangian budget analysis scheme. 
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A. NOGAPS SIMULATIONS 


The "model-data" sets used for this study are from the 
simulations of Sandgathe (1981) to determine the role of 
sensible and latent heat fluxes on extratropical eyed 
genesis. The model integrations started vith idealized 
zonal initial conditions representative of transition season 
climatology in each hemisphere. During the first fev days 
of model integration, a cyclone which resembled a typical 
mid-latitude disturbance developed in each hemisphere. 
After about 6 days of model integration, a variety of secon- 
dary, smaller-scale disturbances developed. Most of these 
secondary cyclones resembled polar lows in that they formed 
and developed in the western (polar) sector of the parent 
low. Some of these disturbances developed rapidly and even- 
tually replaced the parent cyclone as the primary 107 
center. Others migrated cyclonically about the parent Ms 
with little change in intensity. One of the goals ol ma 
diagnostic studies with  numerically-simulated data 154772 
determine quantitatively the mechanisms which cause some of 
these lows to develop whereas others do not. Then the same 
diagnostic analysis will be applied to polar lows chose 
from the FGGE data to demonstrate the validity of conclu- 
sions that are drawn from the model-generated set. 

The NOGAPS model used to simulate the polar lows for 
this study is a state-of-the-art atmospheric model developed 
by the United States Naval Environmental Prediction Reseamep 
Facility (NEPRE), based on a version of the UCLA general 
ClLrcwlearitom model: NOGAPS is a siema-coordinate primi = 
equation model with six vertical layers and rt utm ppc 
Arakawa scheme C staggered grid horizontal differencing 
Scheme. The adiabatic model dynamics are described in detail 
by Arakawa and 1077177777 The model diabatic package 
includes the Arakawa and Schubert (1075 cumulus 


parameterization scheme, and the Randall (1976) version of 
the Deardorff (1972) planetary boundary layer parameteriza- 
tion. A good overview of the model diabatics is given by 
Rosmond (1981); it will not be reproduced here. For the 
experiments of Sandgathe (1981), a sectorized version of 
NOGAPS featuring cyclic east-west continuity was used. A 
Nor zontal domain of 60° longitude by 180° latitude (pole to 
pole) with a resolution of 2.5° longitude by 2° latitude was 
used. The version of the NOGAPS model used by Sandgathe 
(1981) is not the present version. 

Sandgathe (1981) investigated the role of air-sea fluxes 
and latent heating in the development of the secondary 
cyclones by selectively witholding these physical processes 
at various stages of development. In the experiments where 
515500:  surtace moisture flux, or all diabatic processes 
(including latent heating) were witheld, development of the 
secondary lows was hindered. Sandgathe suggested that the 
model development of the secondary lows was brought about by 
low-level thermal processes, and that surface fluxes signif- 
meanitly modify the low-level environment, and Chus 
contribute to development. 

The polar lows in these simulations develop within a 
nanl stic context. Since the model has been integrated for 
6 days before they appear, it is expected that the full 
Spectrum of waves is present. Therefore, the effects of 
nonlinear wave-wave interactions are included, and the polar 
Mes grow in the polar air stream of the parent low, as is 
the case in nature. 

Medel output 15 initially on a staggered grid, and in 
77177 vertical coordinates. Prior to diagnostic calcula- 
tions, the meteorological variables were interpolated onto 
pressure surfaces, smoothed to remove unwanted noise, and 
then interpolated onto the diagnostic budget grid using a 


EE € spline fit to the smoothed isobaric gridpoint data. 
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B. FGGE ANALYSIS 


The First GARP Global Experiment (FGGE) has resulted in 
the most extensive global atmospheric data set ever 
collected by the meteorological community (Halem et “alam 
1056207 Particular emphasis vas placed on data-sparse areas. 
The conventional observational data base comprising surface 
(land and ships) observations, rawindsondes, pilot balla 
and aircraft data vas supplemented by satellite-derived 
soundings and cloud-drift winds, ocean buoys and reconnowm 
sance flights with dropsondes. This improved data base is 
particularly useful for studying maritime weather systems. 

The FGGE data have been combined into dynamically 
consistent analyses at the European Centre for Medium Range 
Weather Forecasts (ECMWF), as described by Bengtsson et al., 
7-7 The ECMWF data-assimilation system consists of a 
three-dimensional multivariate optimum interpolation, a 
nonlinear normal mode initialization and an automatic sus | 
for data checking. A 15-level model with a horizontal reso- 
lution of 1.875? lat/long is used for the dynamical 2S 0 
lation which results in the level III-B analyses used in 
this study. Fields archived are geopotential heights, hori- 
zontal wind components, mean sea-level pressure, initialized 
temperature, relative humidity and initialized vertiz m 
Velocity “tor thew 15 standard levels up to 10 map 
Climatological sea-surface temperature (SST) data are used 
in the physical parameterization package of the model. 
Global analyses for 0000 GMT and 1200 GMT are archived for 
the period December 1978 to 27ü177 77212 During the two 
Special observing periods, January/February and June 2 
additional analyses were produced for 0600 GMT and P9000 M 
Previous studies using the six-hourly FGGE data (Cook 772 
Wash and Cook, 1985; Calland, 1983; Wash and G eee 
noted problems with the data at 0600 GMT and 1800 GMT, 


Therefore, data from these times vere not utilized in the 
present study. 

The FGGE data were interpolated onto the budget grid in 
the same manner as the model-generated data. However, since 
the FGGE data were already on isobaric surfaces, no vertical 


interpolation was required. 


C. QUASI-LAGRANGIAN DIAGNOSTICS 


The quasi-Lagrangian budget approach used in this study 
relates time changes in the volume average of a basic meteo- 
rological quantity (e.g. mass, vorticity, heat, moisture) to 
net gain due to transport of that property through the 
boundaries of the budget volume and to sources and sinks 
within the volume. The budgets presented in the following 
Meeetons are calculated by replacing f in (A.l) (see 
Appendix) by the property of interest. Vertical velocity is 
computed using the kinematic method as described in the 
Appendix. The sources and sinks are identified as dynamical 
and physical processes within the volume, and are determined 
via the governing equations for the property. Thus, the 
budget analysis provides a quantitative measure of cyclone 
development or dissipation in terms of its interaction with 
the environment and in terms of the various physical and 
dynamical processes within the budget volume. 


The budget residual includes the effects of interpola- 


27 and truncation errors, as well as additional physical 
2 “ses not included in the budget formulation. When 
22:57vecd data are used, additional contributions to the 
£E jual arise from inherent data uncertainty, and from 


r Id scale processes not detectable on the discrete grid. 
When model output is used as data, it is possible to formu- 
late the governing equation in a manner that is mathemati- 


cally consistent with model processes. In that case, the 


d S 


budget residual results solely from truncation and interpo- 
lation errors. 

The budget volume is centered on and translates with a 
distinguishable feature of the developing cyclone (or polar 
low). This study utilizes the spherical geometry of the 
quasi-Lagrangian diagnostic budget framework formulated by 
Johnson and Downey (1975). Spherical budget volume geometry 
is especially suitable because many atmospheric systems are 
quasi-circular in shape. The size of the volume is chosen 
ae approximate the scale of the cyclone: The 
quasi-Lagrangian approach permits separation of the advec- 
tive effects associated with the translation of the system 
from the processes which force development. The transport of 
storm properties is determined relative to the moving storm. 
Features of cyclone development are interpreted in a mean- 
ingful manner as changes relative to a moving storm. The 
governing equation for the property of interest is expressed 
in the quasi-Lagrangian framework by replacing the Eulerian , 
local rate of change (23/5t) with S/Gt T 
the rate of change following the system, and € is the trans- 
lation velocity of the system which is computed as descriem 
in the Appendix. 

Budget results are often presented in time-section form. 
Because the budget is evaluated at every model output or 
data assimilation time, these results are used to describe 
the time-evolution of storm processes. A cubic spline =a 
in time is used to provide an instantaneous estimate of the 
time-derivative term. This technique preserves the original 
time resolution of the data by removing the need for time- 
averaging. This feature is especially important for the m 
data which are sampled only every 12 h. Most previous ouc2ə 
studies employ a finite-difference approximation to the eae 
derivative, which is subject to large truncation error SiS 


minimize the time truncation verre: the time difference is 


usually taken over the smallest possible increment and the 
g maining budget terms time-averaged over this period. 

The budget approach emphasizes area and boundary aver- 
aged processes over the budget domain. To examine spatial 
details of storm processes, horizontal fields of terms in 
the governing equation for the property of interest are also 
presented for selected times and levels. These fields are 
also presented in storm-following coordinates. 

A bicubic spline is used to interpolate the data onto 
the budget grid and to compute the various terms for the 
horizontal representation. Because the data are represented 
Böoughout the x, y plane by the bicubic spline function, 
the values on the budget grid are obtained by evaluating the 
255: at the budget gridpoint location. Analytic differen- 
tiation of this cubic yields mutually consistent estimates 
of the horizontal derivatives up to degree two. This is 
computationally efficient because terms involving horizontal 
derivatives occurring in (A.l) or the governing equations 
27 calculated from a single spline fit to the original 
data. inescweonsiderably simplifies computation by avoiding 
the further step of calculating finite differences on the 
budget grid, or alternately computing these on the rectan- 
225577 data grid and interpolating the results onto the 
weet domain. Since the spline is initially fitted to an 
area larger than the budget domain area, there is no degra- 
dation near the boundaries as occurs with finite differ- 
ences. wuse of cübic splines also preserves the 
important smaller-scale synoptic features being studied 
here. For example, мә On involves second 
25 Pivatives of velocity. 1175777 1YEİ st finite “difference 
technique requires three grid points which involves values 
£ sent)ng a distance of 7.5? for a 2.5? grid. uc 
undesirable for the scale of disturbance considered here. 


BE ulating such terms from the original bicubic spline fit 


to the horizontal wind field is computationally efficient 


and results in a much more accurate representations s K |: 


local derivative: Further details are provided in the 
Appendix, vhich includes a derivation Of the 
quasi-Lagrangian vorticity advection term in sphericów 
coordinates. 


The cubic spline has additional mathematical properties 
which make it suitable for use in meteorological analysis 
Because both the gradient and curvature of the bicubic 
spline fit are everywhere continuous, the meteorological 
variable is represented as a smoothly varying property. Tu 
one dimension, the cubic spline results from a minimization 
of the domain integral of the square of the second deriva- 
tive of the property of interest (Marchuk, 1982). For wave- 
like functions, this ensures a preferential attenuation of 
Some analogue of short-wave energy. Such a property MS 
clearly desirable, since the original grid permits only 


wavelengths longer than twice the grid spacing. 


Pil s eCLOGENESTS UNDER STRAIGHT UPPER" FLOW 

HOST s study, the Pacific (or baroclinic) polar low is 
further classified into two types. This classification is 
bused on conditions at polar low formation. The first (Type 
I) appears initially as a shallow baroclinic disturbance 
equatorward of the jet streak under straight or even anticy- 
clonic upper-level flow. The incipient low is associated 
with light surface winds and is located to the upstream of a 
mature cyclone. The low develops in the polar air stream of 
the parent cyclone and increases in both vertical and. hori- 
zontal extent. The polar low develops characteristics of a 
BP s small-scale maritime extratropical cyclone, and 
eventually may replace the primary low as the main circula- 
tion center. The second class of polar low (Type II) is 
distinguished from the first in that it forms on the pole- 
ward side of the upper-level jet ina strong surface polar 
air stream beneath positive vorticity advection (PVA) aloft. 
Examples of Type II disturbances are found in Reed (1979) 
70 “1üllen (1983). Because later stages of development are 
similar for both types, discrimination between types I and 
II is possible only at genesis. 

Over the data-sparse oceans it is not always clear 
whether upper-level PVA is present at early stages of storm 
development. An advantage of model-generated data is that 
the early stages of development can be traced completely to 
determine the relative roles of upper-level and lower-level 


processes. 


Type I polar lows - those developing under straight 
Meer level flow - will be examined in this section. To 
27 L rate the analysis technique, the description of the 
A nostic budget results will be quite detailed. Synoptic 


charts, and mass,  vorticity, heat and moisture 5005465777 
employed in this study. Within this diagnostic framevork, 
the role and importance of latent heat release, and surface 
sensible and moisture fluxes vill be evaluated. An example 
of a Type I polar low has been taken from the Simulated data 
set. The North Atlantic storm analyzed by Cook (1983) and 
Wash and Cook (1985), using the FGGE data, will provi i Eam 
observational example. Cook's mass and vorticity budget 
analysis will be extended to include heat and moisture 


diagnostics: 


A.  SYNOPTIC OVERVIEW 
1. NOGAPS Case 


Geopotential and vorticity maps for selected times 
for the model-generated polar low are shown in Fig. 3.1 ас 
The indicated longitude is based on the 60? cyclic east-west 
continuity used for the NOGAPS simulations. The surface Eon 
originated as an eastward extension from a cyclone situated 
far to the west (Fig. 3.1 a). The central pressure of 
incipient disturbance was around 1009 mb. The initial 277: 
mo ilow res Tk o) is generally anticyclonic and omne 
weakly perturbed. A tongue of weak positive vorticity 
located slightly to the west of the surface distturbamece 
The 700 mb vorticity field for this time (not shown) 1527757 
predominantly undisturbed, which suggests that theo, tr 
bance initially was very shallow. During the first 7590 
(Fig. 3.1 a,b), the surface low becomes embedded in the 
northwesterly flow behind the parent cyclone., mn e 
located to the east, and the central vorticity unen ә 
from 5 x< 10-° “s- tono 310] oe It is at this Stas 
that the system is best classified as a polar low. 11697 
mb pattern for this time (Fig. 3.1 e) indicates cons CEEE 


able amplification of the downstream vorticity mının. 
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Düring the first 24 h, Ер ДОС иси ту minimum 
strengthened from -4 x 10-? s-) to less than -10 x 10-” s-" 
vhile the upstream vorticity maximum increased by less than 
Mabe this amount. This increase in the vorticity variation 
meme the flow results in increasing PVA aloft (Fig. 3.1 
2. 

The polar lov subsequently translated about the 
primary lov and eventually became the main center. At the 
207707 the time period (Fig. 3.1 c), the surface lov is a 
meee, closed circulation with a central pressure of 985 mb. 
Thus, the surface pressure decreased about 18 mb in 30 h. A 
secondary vorticity maximum is evident in this figure at the 
20:00 v At 500 mb (Fig. 3.1 f), the short-wave 
cyclonic feature has expanded to a broadly cyclonic flow. 
Miro rticity maximum at this time is strongly asymmetric, 
with cyclonic extensions to the east and southwest. 

Based on the 300 mb isotach analysis (Fig. 3.2 a-c), 
the polar low formed in the zone between the entrance region 
SO, ms jet streaks. After 24 h, the cyclone is 
located south of a well-developed jet streak. Increasing 
negative shear vorticity downstream from the center (Fig. 


3.1 d, e) is associated with the strengthening of this jet 


Streak. The surface cyclone is located equatorward of the 
1 "or the first 36 h of development. At 276 h, the 
surface cyclone lies directly below the jet axis. At 


пити ју (Fig. 3.2 c), the low is located at the poleward 
5 15 region of a eyclonically curved jet. 

иә nor ion fields (Fig. 3.2 d-£) show that the 
mmm@emeest ascent occurs in the eastern sector of the distur- 
207 and exhibits the characteristic "comma" shape. It is 
noteworthy that almost the entire comma lies equatorward of 
the upper jet axis during early stages of development (Fig. 
EN e). At the end of the period (Fig. 3.2 f), there is a 
secondary weaker ascent maximum to the south of the cyclone 


27 : occurs in the expected location of the cold front. 
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The thermal fields for this case are shown in Fig. 
3.3. The initial surface disturbance (Fig. P K  ин 
to the thermal ridge in a zone of appreciable baroclinity. 
Vith this particular phase relationship, the thermal advec- 
tion patterns tend to amplify both the thermal ridge ahead 
of the lov and the cold trough behind it. A zone of conflu- 
ence ahead of the disturbance results in further strength- 
ening of the horizontal temperature gradient after 24 h 
(Fig. 3.3 b). By the end of the sequence (Fig. 3.3 c), the 
baroclinity is weaker, anda general cooling has occurred, 
especially in the western sector. At 500 mb (Fig. 327770 
the initial temperature gradient in the vicinity of them 
is weak. However, the strong cold advection far to the north 
Serves to enhance the thermal contrast immediately to the 
north of the surface development (Fig. 3.3 e). By the end of 
the sequence, strong cooling has occurred in the west m 
sector (Fig. 2:0 7506 which results in thermal вабре-и- 


Similar to those at the surface at bM ISP 
25 FGGE Case | 


The observational Type I case is the slightly more 
intense North Atlantic storm analyzed by Cook (1963795577 
wash and Cook (1235). A detailed synoptic discussio rioni 
this storm is found in those studies; salient features ormi 
will be presented here. Synoptic charts for this North 
Atlantic disturbance are shown in Figs. 3.4 to Би. 
surface low formed just to the west of Iceland and deepened 
from near 1009 mb to 976 mb as it moved southeastward eOr mm 
the British Isles (Fig. 3.4 a-c). The storm developed nau 
northerly flow around a stationary, occluded system Sorco TE 
off the coast of Norway and eventually replaced this as the 
primary circulation center. The incipient low formed beneath 
a 500 mb ridge (Fig. 3.4 d). However, a short aveo 


that formed about 8° latitude to the northues:5 
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surface disturbance developed concurrently vith the lov- 
level cyclone to become a significant upper trough by the 
end of the period (Fig. 3.4 e,f). 

The initial surface development occurred in a 
thermal ridge to the south of a strong baroclinic Zone (ieee 
ба DE As development progressed, the baroclinic zone 
moved south to a location near the surface disturbance (Fig. 


3.5 b,e). At maturity, very little of the origina epee 


clinity remains (Fig. 3.6 c,f). The surface temperatum 
maxima in Fig. 3.5 a,b occur near the locations Of eee 
Strongest surface pressure gradient. This suggests that 


these locations are warmed by the surface heat fluxes which 
are strongest when the surface wind is strong. In compar- 
ison with the NOGAPS thermal analysis (Fig. 3.3), the FGGE 
surface temperature field is rather noisy. Since the surface 
observations incorporated in the FGGE data set include those 
from a variety of sources including land-based stations, 
buoys and ships, the variability is notio urpi Tines 

The 300 mb jet streak analysis shows that the 
cyclone formed beneath the upper ridge in the right exit 
! jet maximum located to the northwest 
(Fig. 3.6 a). At 0000 GMT 27 January (not shown) thee 


center was located directly below the wind maximum. At later 


region of a 40 m s- 


times (Fig. 3.6 b,c) the storm center was below the Me 
exit jet quadrant. Because the polar low formed beneath the 
right convergent exit region of the upper-level jet, it is 


unlikely that the initial surface development was directly 
related “oO 601-266 :1£ 6 r” 

Satellite imagery during the life of this cy clones 
presented in Cook (1983). This imagery shows that the pola 
low formed in the vicinity of disorganized middle and mne 
level clouds in a generally convective flow. The system 
rapidly acquired an organized multilayered comma-shaped 


configuration accompanied by vigorous deep CconvVect om 
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Following maturation, the demise of the storm was accompa- 
nied by dissipation of this cloud feature. 

The kinematically-derived omega fields for this 
storm (Fig. 3.6 d-f) are rather noisy. The characters 
comma shape and other details of the satellite cloud signa- 
ture are poorly reproduced, which Suggests that the FGGE 
level III-B analyses are only able to describe the coarse 
features of this polar low. At the initial time, ascent 
occurs about 800 km to the north of the disturbance cen 
(Fig. 3.6 d) in conjunction with the left exit region 6279222 
upper-level jet streak in Fig. 3.6 a. Later, a broad area 
of ascent is located to the northeast of the surface low 
(hae. 35 би. 


3) Discussion 


The synoptic features of the model-generated cyclone 
and the FGGE case study of Cook are similar in many 
respects. The essential feature vhich characterizes these 
two cyclones is that both formed under anticyclonic uppoo 
flow on the equatorvard side of an upper wind maximum. 
These synoptic features would appear to place these cyclones 
in the Petterssen et al. (1962) and Petterssen and Smebye 
(1971) type A category. However, upper-level forcing appears 
to play an important role in both cases of cyclogenesrs 000 
both cases, lov-mid tropospheric frontogenetical forcinoldu 
apparent at early stages. In the FGGE example, the mass and 
vorticity diagnostics (presented below) showed that "rau 
circulation increases occurred with the incursion 62705 |. 
forward divergent quadrant of an upper jet streak into the 
budget volume (Cook, 1983). These results suggest WENN 
upper-level forcing may play a greater role in such cases 
than the Petterssen description would imply. 

The detailed diagnostic analysis which follor oee 


highlight many of. these te ü: Similarities between 
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model-generated and FGGE data will demonstrate the realism 
of the simulations and provide further justification for 
analysis of simulated polar lows. Although two very distinct 
storms are studied, the diagnosis of similar physical 
processes vill be regarded as strong support for the impor- 
tance of that process. Differences may be due to dissimilar 
storm processes or may indicate inadequate representation of 


physical processes in the model. 


B. MASS BUDGET ANALYSIS 


The mass budget permits some insight into the strength 
of the vertical mass circulation associated vith these 
developing polar lovs. The mass budget is comprised of net 
horizontal convergence of mass flux vhich is exactly offset 
by vertical flux divergence. These terms are necessarily 
equal because the vertical transport 1S computed from 
kinematically-derived vertical velocity fields as described 
in the Appendix. 

The mass transport is computed for nine 100 mb deep 
İB "rs for a 5” latitude radius (555 km) budget volume. The 
maximum surface absolute vorticity is used to locate the 
budget volume at each time, as described in the Appendix. 
Instantaneous budget terms are evaluated at each data time 
(every 3 h for NOGAPS output and 12 h for FGGE case). 

The time-section of horizontal mass flux convergence for 
the model-generated polar low is shown in Fig. 3.7 a. The 
2” "3€ of 9.69 x 10”” kg is the air mass contained in a 
1 "b deep 5” radius budget cylinder. The lateral transport 
2 2001ts a 2 tvo-layer structure with low-level mass inflow 
2 0v an outflov layer. Meme ee of the convergent 
layer accompanies polar low development with a rise in the 
level of non-divergence (LND) from near 750 mb to 490 mb. 


The rnematically-derived vertical velocities (Fig. 3.7 b) 
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reflect the strength of the area-averaged mass circulation. 
The strongest increases occur during the first 24 h while 
the cyclone is located to the south of the jet, with slow 
strengthening thereafter. 

The mass circulation has an oscillatory fluctuation 
superposed on the slower variations associated with the 
development of this polar low. This cyclic "pattern 
reflected in the surface pressure tendency trace (see A.2 in 
the Appendix) as a 12 h oscillation superposed on the slower 
deepening trend. This oscillation was earlier documented by 
Tallman (1982) who attributed the fluctuations to the pres- 
ence of an external gravity wave propagating within the 
model. Time sections (not shown) of the zonally-averaged 
mass tendency showed that these waves are of meridional 
wavenumber one, and therefore represent a periodic transfer 


of mass between hemispheres. Because this "sloshing" is an 
artifact of the model, it was not analyzed Ii Cen 

For the FGGE case, the mass circulation (EH 777 
develops a two-layer structure similar to that for the 
NÜCAPS” case (pipi 507-001 During the first 12 h orta 
disturbance, the mass circulation is that of an dree ee 
with convergence aloft and divergence near the surface. The 
omega time-section (Fig. 3.8 b) shows that the largest 
increase in the strength of the mass circulation oc w 
between 1200 GMT 26 January and 1200 GMT 27 January MN 
further slower growth until] 0000m CNIT oE nr 7 Some of 
this initial increase inthe vertical mass circula 
occurs while the cyclone is located on the equatorward side 
of the jet streak. Unlike the NOGAPS example where the mass 
circulation remains strong throughout the development 
FGGE case exhibits a decline in the strength of this r A 
lation. This decrease 1S Supported by weakening of the cloud 
features in the satellite imagery toward the end of the 


sequence. 
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25097 ORTECİTY ANALYSIS 


The quasi-Lagrangian isobaric inviscid Moi © re Hy, 
equation 
Š as 
— — E == ә — —— + 
a(S, +f) (V — C) Vic. +f) LU 3p 
A B E 
gu Jw _ ov ow 5: 
üc TY St (35 öy 3p 37) ( ) 


D E 


relates storm-relative changes in absolute vorticity (term 
Pto horizontal and vertical vorticity advection (terms B 
52570) divergence (term D) and twisting (term E) sources. 
This equation is used to specify the source/sink terms 
(terms D and E) for the vorticity budget. The budget equa- 
tion is obtained by integrating each term in (3.l) over the 
pace rec” radiús budget volume which is enclosed between 
pressure levels 1( ) mb apart, “nd partreroning the hori- 
zontal and vertical fluxes into mean and eddy components, as 
described in the Appendix. The budget terms are calculated 
by replacing f in (A.1) by the absolute vorticity, and S by 
HE D and E from (3.1). The terms in (3.1) are also calcu- 


lud directly and displayed on horizontal plots. 
l.  Simulated Case 


Thesevertical time-section of quasi-Lagrangian abso- 
5 vorticity tendency for the NOGXNPS data is shown in Fig. 
7. taste a by, therelis weak “anticyelogenesis aloft. 
Cyretogenesms initially occurs below 600 mb but spreads 
rapidly throughout the depth of the troposphere. While much 
usss increase in the mass circulation occurs during the 
7 24 h, the largest vorticity increase occurs near 300 


£ 7 V . h. These rapid vorticity increases occur as the 


45 


300 mb jet crosses to the south of the surface low center. 
These vorticity changes are described by the various trans- 
port and internal generation processes which are evaluated 


in the vorticity budget described below. 
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bae. ix Time evolution of quasi-Hagranelamn 4550775 


vorticity tendency (10-° s-!) for NOGAPS 
Type I. u 


Vertical profiles of vorticity budget results m 
three stages of polar low development are presented in Fig. 
3.10. A budget term is positive when it represents a process 
which produces a “VOTETCLEY increase. The observed 
quasi-Lagrangian tendency of absolutemvorticity (5052620500 
is equivalent to the sum of vorticity changes due to tchek aks 
terms graphed below, plus a residual. The mean and eddy 
terms cn pcs 3.10 d-i are the mean and eddy flūx con R 
gence terms computed from (A.1). The horizontal eddy termi 
approximately equivalent to the area-averaged advection of 
absolute vorticity. Equality occurs when the convergence and 


mean horizontal vorticity flux convergence terms are equal 
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(see equation in the Appendix). This is evidently the 
case here, except near the surface. Other vorticity budget 
studies in which divergence and vorticity advection parti- 
tions were explicitly computed (e.g. Calland, 1983; Cook, 
8983) yielded lateral eddy flux and vorticity advection 
partitions within a few percent in the mid and upper 
troposphere. 

The existence of small negative residuals indicates 
BENE “cicity Sink not specified in (3.1). Model diffusion 
processes and surface friction probably account for these 
residuals, as suggested by the opposite signs exhibited by 
the relative vorticity and the budget residual (Fig. 3.10 
meee itl anticyclonic areas, diffusion acts as a vorticity 
source by eroding the negative vorticity. The smoothing 
algorithm used prior to interpolation acts as a low-pass 
filter, and would therefore represent an additional diffu- 
sive effect. The large negative residual at final stages of 
mess storm (Fig. 3.10 c) ucun .comıunctlon with large 
relative vorticity, and thus appears to be the result of 
nost diffusion processes. The budget is closed to an 
acceptably small residual, which strengthens conclusions 
regarding the relative importance of the various vorticity 
Memeration or destruction processes. This further reinforces 
275 .ustification for use of model output as data, since 
observational studies using larger budget volumes with time 
averaging have yielded larger residuals (e.g. DiMego and 
25 et, 19982) Cook, 1983, Chen and Bosart, 1979). 

би cüudcet results for the first 21 h of development 
Seeeest that the initial anticyclogenesis aloft is associ- 
ated with upper-level divergence that is offset by a small 
5 “.ive horizontal eddy term (Fig. 3.10 d). Although the 
HD City advection at the initial time is small (see Fig. 
1114), СҺе positive eddy term in Fig. 3.10 d results from 


5 n?”lification of the 500 mb vave that occurs during this 
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21 h period (see Fig. 3.1 e). This growing upper-level wave 


is coupled with the developing surface disturbance. 
Increasing convergence in the lover troposphere gives 
ferticity production there, but this is opposed by a 


negative horizontal eddy term vhich arises from the initial 
asymmetry of the lov. Except near the surface, the net 
contribution from the mean horizontal and vertical flux 


convergence terms (Fig. 3.10 d,g) is negligible. 


At later stages of development, absolute vorticity 
generation becomes significant at all levels. Near the 
surface, vorticity generation occurs as a result of an 


increasing convergence source term and is countered by the 
en C tS of friction, and vertical eddy flux in the volume. 
The horizontal eddy flux term becomes negligible below 700 
mb, which reflects the increasing symmetry of the low-level 
disturbance, and the vertical orientation of the vorticity 
features in this layer. FPriction iS not calculated explic- 
msn this formulation. A frictional sink is inferred from 


B NbUdser residual (Fig. 3.10 b,c), which is predominantly 


negative for the surface layer. There is a net “lossğof 
i evel vorticity by the mean circulation, with vertical 
export exceeding the mean horizontal inflow. However, this 


Mees its not sufficient to cancel the large gain from conver- 
gence. As cyclogenesis occurs, convergence becomes more 


E “nt in producing vorticity growth due to the growing 


Meeolute vorticity coefficient in the convergence term. AS 
the LND increases, vorticity generation by convergence 
тистеазез in vertical extent. j Ser Ph 0.” surface 


vorticity due to low-level convergence was noted in other 
2: . agrangian vorticity budget studies (Chen and Bosart, 
5” bDiMego and Bosart, 1982; Calland, 1983; Cook, 1983: 
£X 7 and Lin, 1984, Vash and Calland, 1985). 

Aloft, the strong cyclogenesis observed near the end 


of the sequence is consistent with the rapid growth of the 


49 


horizontal eddy term (PVA) that is partially offset by the 
divergence sink term (Fig. 3.10 e,f). The tvisting term and 
vertical vorticity transport are important only in em 
middle troposphere. Vorticity generation occurs as a result 
of vertical eddy flux, and is offset by tvisting effects 
(ip. 71005 r E Upvard motion ahead of and subsidence 
behind the budget volume center result in considerable hori- 
zontal gradients of omega, especially near the LND vhere 
omega is largest. Thus, it is in this middle layer that the 
twisting term £ o Eh” DiMego and Bosart (19508 


computed quasi-Lagrangian budgets for the transformation of 


tropical storm Agnes (June, 1972) into an extratropueem 
cyclone. They also found generation by vertical advection 
offset by twisting. A similar relationship between these 


terms was also noted by Calland (1983) and Chen and Bosart 
(1979) for cases of maritime extratropical ecyclogenesis 
which occurred during outbreaks of cold continental arr MD 
would be expected that the twisting term would be a rather 


important vorticity sink in the development of small=s m 


disturbances in a highly sheared environment. Because this 
term depends on the -horizontal gradient of vertida 
velocity, it apparently acts to damp small-secalerbaroc m 


disturbances in the wind field. | 
spatial distributions of the vorti -i generation 
terms in (3.1) for 500 mb are illustrated in PIE OMM н 
quasi-Lagrangian vorticity tendency (g-i of this figur әрк. 
the sum” cof vorticity advection? 32707 convergence (d-f), 
vertical advection and twisting. Because vertical advection 
and twisting effects are smaller in magnitude than the other 


terms and offset one another to a high degree, they are not 
included in Prez pee 


At the initial time, the pattern of vor m 
tendency (Fig. 3.11 g) in the vicinity of theweve tea E 
determined primarily by the convergence term (Fig. 23:11:7022 
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Weak vorticity advection (Fig. 21) reflects the 
initially undisturbed nature of the upper flov and confirms 
this cyclone as Type I. The patterns of divergence (conver- 
gence) occur in conjunction with the upward (donan 
vertical mass circulation east (west) of the disturbance 
center (Fig. 3.2 d), and are largely unrelated to the 
initial vorticity advection pat reue 

The initial pattern of positive vorticity tendency 
to the west and negative to the east results in a vorticity 
pattern that tilts westward with height and leads to the 
development of a vorticity gradient along the upper flow. As 
a result, vorticity advection increases rapidly to become 
the dominant termi(hue. 3 LIIDE At this stage, it 157722 
possible to say whether this cyclone is Type I or 1150162415 
it is developing in conjunction with an amplifying upper- 
level wave. This emphasizes an advantage of model-generated 
data that was noted earlier - namely that we can trace a 
weak disturbance back in time to a point where there is no 
EEE 

The PVA maximum in Fig. 3.11 b is highly corrr 
with the upvard vertical velocity maximum in Fig. 3.2 e. 
Since the thermal gradient at 500 mb (Fig. 3.3 e) at this 
location is weak, it is unlikely that the upvard mor. 
results from warm advection. The amplification of the upper 
waves and the growth of the secondary mass circulation 
appear to be mutually cooperative processes at early stages 
of cyclone development. 

Because vorticity advection occurs dovnstream 22 
the vorticity extrema, subsequent centers of vor r m 
tendency (Fig. 3.11 h) are progressively displaced 46222 
Stream. This tendency eventually leads to an eastward exten- 
sion of the cyclonic vorticity maximum and pushes Cher 


well ahead of the cyclone center (see Fig woe eee 
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By the end of the period, vorticity generation (Fig. 
251 2) occurs in two distinct lobes to the east and south- 
west. The vorticity advection term contributes most to the 
generation. The eastern lobe occurs in conjunction with 
upward motion in that area (Fig. 3.2 f). Convergence at 500 
mb (Fig. 3.11 f) in this region is probably a result of the 
increased depth of the convergent layer at this time (Fig. 
3.7 a) since divergence occurs at 300 mb (not shown). The 
area of positive vorticity tendency to the soùthwest of the 
surface low is the result of upper-level convergence implied 
by the subsidence in that location (see Fig. 3.2 f). This 
tendency results in the southwesterly extension of the 
5 17icity maximum in Fig. 3.1. PVA is very efficient in 
maintaining this lobe, because it occurs near the upper vind 
maximum. The area of weak ascent ahead of this region 
results in an area of low-level convergence to the south of 
the surface low center. This convergence is responsible for 
2” ГогпабСјоп of the secondary vorticity disturbance "w" in 
2- :.1 c. 


7 FGGE Case 


mh elat ive vorticity tendency (Fig. 3.12) exhibits 
large increases throughout a deep layer centered on 0600 GMT 
57 Panuary, ufenhodocreasino” vorticity after 1200 GMT 28 
Enuary. The strongest  cyclogenesis occurred during the 
period when the surface cyclone crossed to the poleward side 
of the jet axis. 

Vorticity budgets for the FGGE case are presented in 
IB 3.13. Lunuussuno first 12 nh, the calculated budget 
2— (o not account for the lov-mid tropospheric cyclo- 
genesis. The budget terms at this early stage indicate 
convergence aloft above lov-level divergence, vhich are 
£ "stent vith the subsidence during this period (Fig. 3.8 


b. Bu iirec budget residual indicates that an additional 
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Fig. 3.12. As in Fig.-3.9 except or TOCHTE CE 


vorticity source is needed to explain the positive vort Icas 
tendencies which occur in conjunction with a vorticity e ии 
due to divergence and eddy transport losses. This residual 
possibly arises from data inconsistencies. During the 11572 
12 h of development, the elevated terrain of Greenland is 
included in the western part of the budget area. In winter, 
radiational cooling results in a cold layer near” 07 
surface. The surface vinds in this layer are strongly decou- 
pled from the synoptic pressure patterns and are represe 
tive of a very shallov surface layer. Over elevated terrain. 
the ECMVF analysis extrapolates to pressure levels belov the 
topography. Thus, the terrain-dependent Surface wind infor- 
mation may be applied to a much deeper layer than is appli- 
cable, and this may account for the large budget residua 3 


at least bêlow about UU ии 
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Between 0000 GMT January 27 and 0000 GMT January 28, 
convergence near the surface causes a  vorticity increase. 
This is offset by friction (negative residual) in the lowest 
layer. Generation aloft is from PVA and is partly offset by 
divergence. During the final 24 h of the sequence, 
decreasing  vorticity occurs throughout the whole  tropo- 
Sphere. Aloft, this is largely a result of a much reduced 
eddy flux term which is caused by the decreasing vertical 
tilt of the system at this time (see Fig. 3. c; f) eR 
the surface, decreasing low-level mass inflov inhibits 
further generation. As for the NOGAPS case, the horizontal 
and vertical mean flux convergence terms (Figs. 3.13 d-f and 
g-i) are equal and opposite so the net effect of these terms 
is negligible. 

Horizontal depiction of the terms in (3.1) at 507902 
is given in Fig. 3.14. Initially, weak NVÀ occurs direc 
above the surface low (Fig. 3.14 a), which shows that sa 
is a Type I polar low Thot initiated by PVA. The area of PVA 
to the east is associated with the vorticity minimum near 
that Location Seco ə 0 s Generation by convergence 
occurs above and to the northeast of the low center Gs 
3.14 d). The area of strong divergence and PVA to the moisi 
of the center is associated with the left exit region of the 
jet streak to the northwest (Fip 70067 Since the jet- 
related divergence is located about 800 km from the surface 
center, it is unlikely to be responsible for genesis of this 
Storm. 

At 1200 GMT 27 January, the strongest positive 
vorticity tendency occurs to the south and vest 0770202 
center (Fig. 3.14 h), with negative tendencies to the они 
and east. This is consistent with both the southward mays 
ment Of the upper Wer streak. and the growing amplitude of 
the 500 mb trough/ridge in Fig. 3.5 e. As the cyclone Tames 


occludes, the largest positive tendency (Fig. 3.14 $9906 
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in lobes to the east and the south of the lov center (note 
that the center of Fig. 3.14 at 1200 GMT 28 тбабт әри. 
located vell to the vest of the surface lov center in Fig. 
3249). This is similar to the final vorticity tende mieki 
pattern for the NOGAPS case in Erg ири 


3. Summary and Discussion of Vorticity ms 


The development of these two Type I cyclones has 
been described in terms of the various vorticity generation 
mechanisms. Mass circulation increases lead to von enSa 
generation from the divergence term. The resulting grovth of 
upper-level perturbations causes increasing vorticity advec- 
tion in conjunction with the divergent secondary circula- 
tion. Vorticity advection leads to deformation of the upper 
vorticity patterns. In both cases, strongest cyclogen: ди. 
aloft appears to immediately follow the time when the 300 mb 
jet streak first crosses to the south of the low cen a 
Surface vorticity generation by convergence maintains the 
surface cyclone against the effects of frictional diss 
tion. The nonlinear deformation of the upper-level 
vorticity field appears to create an environment favorable 
for the next secondary cyclogenesis. 

As the low develops, the growing vertical motion 
results in increasing upper-level divergence. From (3.1) 
this may be written 

кә ии (3.2) 


(EIE 


which neglects the twisting and vertical advection (2008 
that were shown in the vorticity budget analysıs tol: 
“malı By definition Type I polar lows form under €5502 
tially straight upper-level flow so the initiai D 
small. From (3.2), the growth of upper-level divergence ls 


therefore balanced by decreasing vorticity. Since the 
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Mtia] location of the polar low is in the lov absolute 
vorticity environment of the warm side of the jet streak, 
5 reduced denominator of (3.2) means that only slight 
anticyclogenesis is required to balance the mass circulation 
increases. Despite growth of upper-level divergence, very 
r le PVA is initially required for vorticity balance. 

The FGGE and the model-generated polar lows exhibit 
some differences. In the FGGE case, Vorticity generation 
2” rs throughout a deep layer during the first 24 h (Fig. 
5:2). Maximum cyclogenesis coincides with tthe end of the 
Strongest mass circulation increase, and decreases there- 
after. In the NOGAPS case, maximum cyclogenesis aloft occurs 
well after the initial mass circulation increase. One reason 
for the later cyclogenesis in the NOGAPS case is the intense 
downstream ridging aloft which offsets vorticity increases 
in the trough early in the development. Another reason is 
the longer period that the NOGAPS case remains in the lower 
vorticity environment equatorward of the upper-level jet. 

The decrease of the mass circulation in the later 
half of the FGGE case gives rise to weakening low-level 
convergence and decreasing vorticity generation. In the 
NOGAPS case, the maintenance of the vertical mass circula- 
tion ensures continued vorticity generation from the nonli- 
near convergence source term. TE will remain to De 
determined why the strong mass circulation persists longer 
in the model-generated case. 

For the FGGE case, Cook (1983) shoved that the most 
B5- oTrous vorticity increases were associated with the incur- 
255 0 0£ the forvard divergence quadrant of a jet streak into 
the budget volume (see Fig. 57700: Initial development of 
27”  .sturbance cannot be attributed to jet-related mass 
Mieeeulation increases, since the low formed equatorward of 
Meee jet streak (Fig. 3.6 a) under an area of weak NVA aloft. 


The upper-level jet streak was initially uncoupled to the 


5- 


low-level disturbance. However, once the forward divergence 
quadrant of the jet became superposed vith the surface 
disturbance, Strong coupling and vigorous cyclogenesis 
occurred. 

It is important to realize that the vorticity 50 E 
results do not imply a mechanism for cyclogenesis. Rather, 
vorticity changes are viewed in a balance context. The 
actual changes in vorticity reflect the continual hnydr2 
tatic adjustment of geopotential to changes in the thermal 
structure, vhich are themselves dependent on a variety of 
physical processes. Thus, the vorticity budget results alone 
are unable to yield a complete description of the cyclone 
development. For this reason, a diagnostic analysis of the 


thermal processes was also performed. 


D. THERMODYNAMIC ANALYSIS 


1. Heat Budget Analysis 


The quasi-Lagrangian budget of potential temperature 
relates changes in the volume-averaged potential temperature 


to the various heat flux convergence terms evaluated using 


55: and to the diabatic source and Sink 7 The 
diabatic terms arise from air-sea heat glu radiative 
fluxes and cCondensational heating. In terms of potent magi 


temperature, O9 the thermodynamic eguak momen 


25 EE F. : 
— = -(V-C)'VB8 - u — + (1000, p o 


) (352 

Pu jp E 

s p 
where Q/Cp is the total diabati€ 1641507 vhich is 
specified as a model output yar r yb e P Therefore, it p 


possible to specify completely the model heat budget, with 
budget residuals representing truncation and interpo See 


errors only. 
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Vertical profiles of the heat budget terms for the 
NOGAPS case are graphed in Fig. 3.15. The calculations are 
applied to the same budget volumes as the mass and vorticity 
budgets. The budget is computed by replacing f in (A.1) by 
potential temperature. The diabatic source term is specified 
Merche last term in (3.3). The quasi-Lagrangian potential 
temperature tendency (Fig. 3.15 a-c) is equal to the sum of 
the net mean and eddy flux convergence terms, and the total 
EGbatic heating (Fig. 3.15 d-£), plus the budget residual. 
The horizontal and vertical components of the net eddy 
E£ “tition, mid them cumulus "heatene (part of the total 
heating) are included for convenience as g-i of this figure. 
The horizontal and vertical components of the net mean heat 
flux convergence are each an order of magnitude larger 
(100K /day) than the graphed terms, and are approximately 
equal and opposite. Therefore, only their sum is presented. 

The residual term is generally small relative to the 
leading terms. However, large changes in the vertical 
gradient of potential əmperature near the tropopause cause 
w em ic errors in tie vertical interpolation of tempera- 
ture from sigma to pressure surfaces. The erroneous theta 
values lead to budget residuals that 7 highly correlated 
vith the large flux convergence term (compare the residual 
and the net mean flux convergence profiles in Fig. 3.15). 
Negligible error is introduced from interpolation of the 
momentum terms because the mass circulation is balanced to 
221 precision at each gridpoint. 

At early stages of polar low development, potential 
temperature changes are slight (Fig. 3.15 a). However, as 
puesevyclone grows, net cooling occurs (Fig. 3.15 b,c) as the 
low moves toward higher latitudes. Tufs?ool ne 15 Consis- 
tent vith the rapid cyclogenesis aloft later in the develop- 
EU Fig. 3.9). The largest terms in the heat budget are 


diabatic heating and the approximately equal and opposite 
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adiabatic cooling that is associated with the mean mass 
Es lation Fig 3.15 d-f). Diabatic heating arises prima- 
rily from condensation and latent heat release in ascending 
dır. 

Compared vith the tvo leading terms, the potential 
temperature tendencies, and horizontal and vertical advec- 
tion (eddy terms) are small in this area-averaged setting, 
except near the tropopause. The horizontal eddy term is 
small, since the 5° latitude radius budget area includes 
both strong cold and warm advection. To relate thermodynamic 
terms to synoptic features of this developing storm, heat 
budget terms are depicted in a horizontal format in the next 


“ubsection. 


2. Thickness Tendency Analysis 


Maintenance of hydrostatic conditions is an impor- 
tant physical constraint in the atmosphere, and is guaran- 
teed in the numerical model. Temperature increases 


(decreases) in an atmospheric column must be accompanied by 


increasing (decreasing) geopotential thickness for the 
layer. For example, net temperature decreases in the column 
above the developing polar low observed in Fig. 3.15 c are 


compatible with the vorticity tendency profile for the same 
0” oc (Fig. 3.10 c), which indicates stronger cyclo- 
genesis aloft compared to that near the surface. 

To gain more insight into the relative importance of 
the individual thermodynamic processes in polar low develop- 
Ep the terms in the quasi-Lagrangian thickness tendency 


equation have been evaluated as follows: 





E. 900 из 
тә = | 1-(V-€)-VT e (Q--52)w «D ]d(in p) (3.4) 
i 500 p P D 

A B C D 


63 


The layer below 500 mb vas chosen for this calculation 
because it is most relevant for polar low development. Also, 
interpolation errors in the temperature values are larger 
aloft. An average potential temperature increase of 1°K for 
this layer corresponds to a thickness increase of approxi- 
mately 15 gpm. The quasi-Lagrangian correction to the 
thermal advection (term B) is employed. The vertical advec- 
tion (term C) is proportional to the product of?” vere 
velocity and static stability: This term also will be 
referred to as adiabatic cooling in ascent regions (warming 
in descent areas). The diabatic heating (term D) is the 
vertical integral of the model-derived total heating. The 
net quasi-Lagrangian thickness change (term A) is calculated 
as the sum of the three terms in (3.4). Reasonable agree- 
ment was found between this sum and the explicitly calcu- 
lated thickness tendency field (not sizim E The latter 
involves evaluating time derivatives at each gridpoint. It 
was concluded that the sum of the three terms in (3.4) 
yields reliable estimates of the thickness tendency. 

The thickness tendency equation analysis for the 
NOGAPS case extends a similar analysis by Sandgathe (1981) 
obanı Clow nis 7 His analysis used the Eulerian 
form of the thermodynamic equation in sigma coordinat es AS 
present technique employs the quasi-Lagrangian adju.tment to 
the thermal advection term, and is presented in isobaric 
coordinates for three stages of cyclone development (Faces 
5057 

During polar low genesisği the thermal advection 
patterns (Fig. 3.16 a) have a large cooling tendency tome 
North: of -Ehe suntace cenuer- This cooling is partiy or oE 
by subsidence varming (Fig. 3.16 d). Diabatic heating (Fig. 
3.16 g) and adiabatic cooling (Fig. 3.16 d) almost rec 
cancel to the east and south of the development. The net 


effect (Pip GE NEED is a tendency for the thin 
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gradient to increase in the area to the north of the surface 
cyclone. This frontogenetic tendency was previously noted 
in the synoptic discussion (see Fig. 3.3). Vertical asen 
ahead of the surface disturbance appears to be initiated by 


weak warm advection and augmented by the latent heat 


release, since the fields of ascent, warm advection and 
diabatic heating are approximately co-located, and the net 
temperature tendency is comparatively small. It is Cimi 


growing vertical circulation which is responsible 7017527 
vorticity generation discussed earlier. 

After 24 h, a strong dipolar pattern of advective 
cooling to the rear and warming ahead of the lov has devel- 
oped (Fig. 3.16 b). In the ascent region, the net effecmuam 
diabatic heating (Fig. 3.16 h) and adiabatic cooling SR 
3.16 e) is slicht C oping, so net thickness changes (Fig. 
3.16 k) approximately reflect the contribution from th€21224 
advection. 

At the end of the sequence, the net warming patterns 
(Fis: SEM reflect a cyclonic rotation of ear MEEN 
patterns. This rotation reflects the deformation Of iia 
thermal fields by the advection term. At this .time? 777 
trend is to eliminate the north-south temperature and geopo- 
tential gradient near the location of the polar lo pp K 
also has the effect of producing a more vertical geopomena 
tial structure. On the south side of the area of coo mame 
the lower left of Fig. 3.16 1, the trend is frontogene miksi 
This trend, and the frontolytic trend near the low cem m 
are consistent with the movement of the jet streak to the 
south of the surface low center (Fig qM m 

Very good qualitative agreement exists between the 
present analysis at 264 h and that of Sandgathe (1931777006 
main difference is the presence of an area of strong mm 
advection to the northwest of the surface eyclone” 115577 
quasi-Lagrangian analysis (Fig. 3.16 b). This difference 


reflects the effect of storm tr u s us SAR 
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2750700 5 required. in the interpretation of the 
thickness tendencies in Fig. 3.16, as geopotential changes 
aloft and near the surface are included. In general, poten- 
tial temperature tendencies are larger aloft than near the 
surface (Fig. 3.15 a-c), as Specific volume increases with 
height. It follows from (3.4) that geopotential changes at 
500 mb account for most of the diagnosed thickness varia- 
tions. Exceptions to this may occur in regions where net 
lov-level heating (cooling) is much larger in magnitude than 
meee aloft. In such cases, most of the thickness increase 
(decrease) would be reflected by falling (rising) surface 
geopotential. 

Comparing Figs. 3.16 k and 3.1 e, it 1S apparent 
that thickness increases (decreases) occur near the 500 mb 
206 (trough). Such localized changes in the thermal 
structure of the troposphere tend to amplify the 500 mb 
emeh- ridge couplet, which is compatible with the progres- 
sive increase in the variation of vorticity along the flow 
noted earlier (see Figs. 3.1d and e). The increased PVA 
implied by this change requires an increase in the upper- 
level divergence betveen the trough and the ridge. The 
associated upward motion cools the layer, offsets the warm 
advection and opposes further amplification of the upper 
wave. However, release of latent heat in the ascending air 
largely offsets this limiting effect. Warming beneath the 
dee continues and leads to further amplification of the 
upper wave. The progressive growth of upper-level divergence 
favors increasing low-level convergence, vurcn leads -to 
məsonsification of the surface cyclone. io pi 1 ЕДис e 
Eun turn strengthens the thermal advection fields and 
" 4165 in further amplification of the upper waves. TRGS., 


1 


Ehe so-called "self-development" Mechanism Tor cyclone 
22 vification proposed by Sutcliffe and Forsdyke (1950) 
Eur: to apply for maritime cyclogenesis on this scale of 


Motion as well. 
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At the -initial Eime” the net diagnosed thickness 
tendency (Fig. 3.16 j) is not conducive to self-development. 
Although the thermal advection at this time produces a weak 
dipolar pattern of the right phase for self-development, it 
is more than offset by the effects of adiabatic cooling. 
Similar results were obtained when (3.4) was evaluated for 
the layer between 900 and 800 mb. Amplification Of (james 
thermal advection dipole caused by development of the low- 
level cyclone eventually results in net thickness changes of 
the correct phase for self-development (Fig. 3.16 k). The 
Strengthening of the thermal gradient to the north of the 
cyclone may also contribute to growth of the thermal advec- 
tion term: 

Self-development is favored by some combination of 
reduced static stability, high relative Tüm EE n and 
enhanced baroclinity. In very stable air, adiabatic cooling 
(term C of (3.4)) is large. Therefore, for self-development 
EG, CeCuc. larger contributions from thermal advection wae 


latent heat release are necessary to offset the limiting 


effect. This is achieved when the thermal gradient is large 
and the air is moist. In a lower static stability envia 
ment, the limiting effect is smaller, and self-development 


will occur for smaller thermal gradients andam r 
CONTENE: 

The Richardson number measures the ratio of static 
stabilty to baroclinity. Therefore, a low Richardson mun 
is likely to be conducive to self-development. Following 


Gall (1976a), a spherical Richardson number is defined» 
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where f,, is the value of f at 45° N. The term involving 
tan@ in the denominator of (3.5) normalizes the vertical 
wind shear so that a given shear is equivalent to the same 
horizontal temperature gradient, irrespective of latitude. 
This formulation is similar to the zonal Richardson number 
25 Call (1976a). According to Gall (1976a) and Staley and 
2701 (1977), reduced Ri favors development of short-vave 
baroclinic disturbances. The reduced values of 700 mb Ri 
(Fig. 3.17) which occur slightly downstream of the distur- 
bance center are consistent with the hypothesis that these 
polar lows grow through baroclinic processes. Additional 
areas of reduced Ri occur to the east and southwest of the 


Mew center (Fig. 3.17 c). 





75 .-.1/7. 700 mb EE (contours every 60 dams EN 
and Richardson number for NOGAPS Type I. Values 
above 10 are not contoured. 


The rotation of the thickness change dipoles is 
Important for later secondary cyclogenesis. Following devel- 


opment of the cyclone, cooling occurs in the southwest 
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quadrant (Fig. 3.161). This cooling, which IS e MI 
result of thermal advection (Fig. 3.16 e758 is consistent 
with the cyclonic vorticity generation at 500 mb (Fig 3 
ı)” occurrıng in Cha: on Furthermore, the patterns of 
cooling at the location "w" in Fig. 3.16 1 are fronte мә 
etic, which results in an area of enhanced baroclinity near 


t 11 


w in Fig. 3.3 c. Thus, the secondary vorticity maximums 
"w" in Fig. 3.1 c appears to be the result of baroc miN 
processes, as Suggested by the Ri minimum near that location 


(ETE... cS 


3. Influence of Latent Heat Release 


eS  —— FT 


Condensation in the ascending air ahead of the 
developing cyclone results in release of latent heat. In the 
thickness tendency analysis presented earlier, diabatic 
heating was shown to promote continued self-development of 
the polar low by opposing the limiting effect of adiabatic 
Goo lamer Fe. sl D. This means that vertical ascent can 
increase further than would occur if diabatic heating were 
absent. This augmented vertical velocity leads to increased 


generation of low-level vorticity and results ina more 


vigorous cyclone. The enhancement of ascent by latent 
heating has been noted by many investigators (e.g, 
Smagorinsky, 1956; Aubert, 1957;  Danarsd we ии Dana ni 


(1964) solved the diagnostic quasi-geostrophic omega equa- 
tion with and without latent heat release for a major m K 
cyclone over the central United States. He found an appr = 


mate doubling or the vertical ascent when latent heating was 


included, which led to strengthening of the low-level 
convergence and high-level divergence. This resulted in an 
increased production of low-level vorCicityE Associated 


changes in the ageostrophic (divergent) wind components 
contributed to kinetic enerey (KE) generation at both Him 


and lov levels. 
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Duma role tor the internal diabatic heat source 
is generation of available potential energy (APE), which may 
be converted to storm KE. The concept of APE was introduced 
by Lorenz (1955), who defined it as the difference between 
the total potential energy of the atmosphere at a given 
time, and a hypothetical state obtained by adiabatically 
rearranging the thermal structure to a statically stable 
barotropic state. As such, it represents the part of the 
potential energy which is available for conversion to KE. 
The existence of APE generation is not a guarantee of an 
increase of storm KE; rather, it indicates an energy source 
potentially available for conversion to KE. Conversion is 
indicated by a positive covariance between specific volume 
m@a@evertical motion (warm air rising, cold air sinking). As 
a first approximation, generation of eddy APE depends on 
covariances of heating and temperature deviations on an 
isobaric surface. 2170277070 0 969), -and Bullock and 
Eobhnson (1971) used the exact expression for APE to show 
that a combination of heating at low-mid tropospheric levels 
in the warm sector, and cooling at mid-high levels in the 
cold sector is most conducive to eddy APE generation. 

Production of APE for the model-generated polar low 
is evident in the previous analyses. The largest diabatic 
Meeting (Fig. 3.16 g-i) occurs above the location of maximum 
positive low-level temperature perturbation (Fig. 3.3 a-c). 
The thermal patterns at 700 and 850 mb (not shown) re 
BDDroximately in phase with those at the surface. As self- 
development progresses, the net effect (Fig. 3.16 k-1) of 
E stic heating, adiabatic cooling and thermal advection is 
to amplify the thermal vave. Moreover, the diabatic heating 
mü ucreasing vith time. The correlation of the diabatic 
heating with the growing temperature wave ensures increasing 


— nration of APE as the low intensifies. 
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Covariances of omega and specific volume evaluated 
at 100 mb intervals have been used to diagnose baroclinic 
conversion (Fiz.777777 Initially, weak eddy APE to KE 
transformations occur in a shallow layer below 750 mb, with 
another layer of weak positive conversion aloft. During the 
first 12 h, the shallow surface layer of positive conversion 
extends to 500 mb and strengthens. A thermally indirect 
circulation (negative conversion) develops above this level. 
Later, positive conversion occurs throughout the wr 
Sphere. The time section of diabatic heating (Fig. 3.18 a) 
is rather similar to that of vertical velocity (Fig 777770 
This suggests that vertical ascent is enhanced by latent 
heat release. The similarity between the ascent and the 
diabatic heating patterns suggests that the temporal varia- 
tion of the APE generation may resemble Fig. 3.18 b. 

This analysis has shown that diabatic heating due to 
condensation in the warm sector enhances both generation of 
eddy APE and conversion to KE in the polar low. This 1522902 
accord with previous work for mid-latitude cyclones ( ири 
Danard) 17:97 Bullock and Johnson. mE Tn Gall, 1928 
Vincent et al., EO T: Robertson and msn s n 1983298 
Furthermore, it is apparent that net warming patterns favor- 
able for self-development (warming in warm sector, cooling 
in cold) imply increasing APE generation and conversions 
Storm kE: Hodur (1934) compared adiabatic and diabatic 
Simulations of a wavenumber 20 (wavelength about 1400 km) 
barocliınıc dissurbumnee” The disturbance grev from a 20222 
initial state typical of conditions during the "senes mM 
Pacific polar lows (Mullen, 1979). Hodur found that ita MEME 
sion of diabatic effects increased the KE of the disturb ии 
by an order of magnitude, as rapid intensification resum 
from a much strengthened baroclinic conversion of eda m m 
to KE. The depth of the layer of positive eddy conversum 
was much greater for the moist experiment which resulted 


a deeper disturbance. 
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The Pacific type polar low is generally a deep 
disturbance having Significant geopotential amplitude near 
the tropopause (Reed, 1979; Mullen, 1979; Locatelli et mm 
1982). At genesis, . the Type I polar low is shallow, Мари 
forms under straight upper-level flow. The release of 
latent heat plays a crucial role in the formation of UPON 
level geopotential waves by contributing to the increase in 
the temperature perturbation in the warm sector at mid- 
levels. The maximum eddy conversion is shifted upward which 
results in a deeper disturbance than would be the case in 
the absence of latent heating (Gall, 1976b). Alternately, 
diabatic heating allows self-amplification of the upper 
waves to progress further. 

Tractons Gio o) suggested that the latent heat 
release associated with precipitation creates the imbalance 
necessary for the onset of development prior to the time 
when vorticity advection alone would initiate development. 
This is certainly the case here, since the initial EY. 
minimal. In this case, the development of upper-level PVA 
can be viewed as the hydrostatic response of the system to 
thermodynamic fererne and ls (nett a Precurs om QUE 


cyclogenesis. 
&. Moisture Budget Analysis 


In the previous subsection, it was shown thia e 
thermodynamic balance of the developing system depends 
crucially on the strength of the latent heat relea s Di 
examine the processes which ensure the supply of water vapor 
to the ascending air, a quasi-Lagrangian moisture budget 
been performed (Fig. 3.19) for the NOGAPS-simulated polar 
low. The moisture conservation equation 

ј 


= = ER E - +T ee p ) (3402 
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expresses the difference between evaporation (E) and precip- 


Econ (P) as the sum of horizontal and vertical flux 


divergences of specific humidity (q) and the 
quasi-Lagrangian tendency term (change in storage). EZ) 
computed as a budget residual. Large accumulations of vater 


vapor result from horizontal and vertical flux convergence 


(Fig. 3.19 d-f). Because there is little change in the mean 


vapor content of the column, there must be a net loss of 
Water from the system (precipitation), as indicated by the 
negative residual for the volume. It should be remembered 


that evaporation of the condensed water is also occurring at 
some gridpoints within the volume. The negative residual 
indicates that precipitation exceeds evaporation when aver- 
aged over the budget volume. 

loss rna water Vapor accumulation occurs as a 
5 JC of the mean circulation. Horizontal moisture flux 


convergence is strong near the surface and decreases to near 


zero above 700 mb (Fig. 3.19 d-f). Near the surfac: o this 
| 

fD offset by vertical tlux divergence. In the mid ar ` upper 

troposphere, most of :he moisture accumulation for precipi- 


tation arises from vertical flux convergence. 
Be Role of Air-Sea Fluxes 


The boundary fluxes of moisture and heat arise from 
air-sea moisture and temperature differences in the presence 
boundary layer winds. In the NOGAPS simulations, the 
model sea-surface temperature (SST) is independent of longi- 
tude vith the maximum SST gradients betveen 45 and 50° N and 
between 32 and 42” S. 

In NOGAPS, the PBL fluxes are assumed to heat the 
lowest sigma layer of the model, which is about 200 mb deep. 
27  rface heat flux of 100 W m-^ is equivalent to a heating 
Hate ot about 4°K per day for a 200 mb deep layer. The 


E sea sensible heat flux is shown in Fig. 3.20. At early 
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stages of development when the surface winds are light, the 


mixes are negligible (Fig. 3.20 a). As the disturbance 
becomes more intense, the fluxes (Fig. 3.20 b) tend to warm 
2:7” coldest air and cool the warmest air (compare Fig. 3.3 


b). These sensible heating patterns oppose the effects of 
thermal advection and inhibit grovth of the thermal vave. 
This tends to destroy eddy APE and oppose self-development. 
Hovever, even the strongest flux in Fig. 3.20 amounts to 
only a few degrees per day of heating. A surface flux of 100 
W m-? causes a thickness increase of only 25 gpm per day in 
75— 0:0 mb surface layer, which is only a small fraction of 
the diagnosed thickness tendency components in Fig. 3.16. 
Thus, in this case, the surface heat flux offers minimal 


opposition to cyclogenesis. 
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6. YFGGE Case 


The budgets of heat and moisture, and the thickness 
tendency analysis have been also performed for the FGGE 
pölarilow. The heat budget results are presented in Fig. 
3.21. The diabatic heating is inferred as the budget resi- 
dual, which also includes computational uncertainty. At the 
initlal time, there is a deep layer of cooling above a 
shallov surface layer of heating (Fig. 3.21 39m It e 
unlikely that a deep layer of diabatic cooling ot "Siem 
intensity would result from radiation or evaporat mog 
Similarly, the layer of diabatic heating of 10° per rm 
the stratosphere at this time is also unrealistic. Thus vine 
diagnosis of heat budget terms using FGGE data is subject 
much uncertainty. It is possible that strong radiation 
cooling which exists near the ice-covered, elevated terrain 
of Greenland during the long polar night may contribute to 
the negative residual at mid-levels. However, the shallow 
layer of sürface heatine any Pie - 3.21 a appears to contem 
dict this explanation. This layer of heating 1s probably game 
to increasing surface fluxes as the disturbance moves || 
ice-covered terrain over warmer water and intensifies. 

The inability to close the FGGE heat budget to an 
acceptably small residual appears to be a result ора. 
uncertainty. Additional error may be introduced by them 
of the archived temperature data which are from the ECMWF 
forecast model first guess (Paegle, 1983). These tempera- 
tures are not in hydrostatic balance with thema s || 
height fields. Wash and Calland (1985) compared |) E 
temperatures with those computed hydrostatically from the 
height fields. They found RMS differences in the tropo 2 
of less than 2°C between the two estimates, while e ma 
Stratosphere, these differences were nearly + CIS 


departure in the stratosphere possibly accounts ro 
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large variability above 300 mb in the residual term in Figs. 
9.21 amc 

At later stages of development, diabatic heating 
(inferred from the heat budget residual) occurs throughout a 


deep layer (Fig. 3.21 b) and is countered by adiabatic 


cooling from the net mean fluxes İ(Fig: 721777 This 
similar to the heat budget result for the NOGAPS case (Fig. 
25215)” In the FGGE case, the strength of the diabatic 


heating near the surface indicates large surface sensible 
heat flux. 

The profiles of temperature tendency during the 
first O ao 2515006) are consistent vith the rapıd 
decrease in static stability during cyclogenesis which was 
discussed by Cook (19837. Folloving Sandgathe (1981), Cook 
used a stability index obtained by subtracting the 1000 mb 
potential temperature from that at 500 mb and dividing Soi 
the pressure difference. For the 4° radius budget volume 
used by Cook, this index decreased from 4.8°K per 100 mb at 
1200 GMT 26 January to a minimum of 2.1?*K per 100 mb a 
GMT 29 January. This minimum approaches the moist adiab sm 
value based on the area-averaged surface temperature at this 
time. For the NOGAPS case, the stability index is appro Si 
mately 4.7°K per 100 mb. This is close to the moist adas 
batic value OLA g k per 100 mb. derived from the 
area-averaged surface temperature of 15°C. This isi 2 
tent with the near-cancellation of the diabatic heatınEa 22 
the adırıbatıc cool IE 

The time section of diabatic heating (Fig. 2 
shows a deepening layer of heating with larger values nearer 
the surface. Strong surface heating is due to heat flu |: 
from the sea. At the initial time, the surface har mi G 
is weak, and the fluxes are small, in spite of large in m 
air-sea temperature differences (Cook, 1983). However, “as 


the cyclone develops, these fluxes rapidly become large and 
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result in an excess of 8? per day of surface heating at 1200 
GMT 28 Qanuary. The later decrease of surface sensible 
heating is due to movement of the storm over Northern 
Ireland, with a large fraction of the budget volume over 
land. The SST gradient also is weaker’ south of 60? N. 

Another factor which may account for the low eleva- 
tion of the strongest heating is the initial vertical 
distribution of moisture. Because the air originated over 
ice-covered terrain, it would be expected to be rather dry 
emor to polar low formation. Large surface fluxes of heat 
and moisture quickly develop as this cold, dry air streams 
southward over the open water. Therefore, most of the mois- 
ture at early stages of development is confined to a shallow 
layer near the surface. Because of this moisture distribu- 
tion, most of the latent heat release associated with ascent 
occurs at low elevations. This deepening layer of low-level 
heating contributes to the rapid destabilization described 
earlier. 

ME ro O conyersion term (Fig. 3.22 b) is 
mostly positive above 800 mb, which supports the hypothesis 
Meare this cyclone grows through baroclinic processes. 
However, this term is smaller than for the NOGAPS case, and 
there is a layer of negative conversion near the surface. 
This reduced conversion is probably due to the much larger 
Surface heat flux in the FGGE case, since strong heating in 
the cold sector tends to reverse the phase of the tempera- 
ture wave near the surface. Thus, It appears that the 
2 (— sensible heat flux directly opposes baroclinic 
development. 

The FGGE budget of vater vapor is presented in Fig. 
5. At early stages of development, a deep layer of 
significant evaporation is implied by the positive budget 
residual. such evaporation is unlikely except near the 


surface where cold, dry air moves over progressively warmer 


Ə 1 


HEAT BUDGET RESIDUAL (°C per day) 


150 
290 
350 
+50 
550 


650 


Pressure (mb) 


850 





950 
26 20 28 29 


TIME (January 1979) 


Pet 


Pressure (mb) 





26 2 28 29 
TIME (January 1979) 


bığ” ај Heat budget residual for EGGE Type 1. 
b) As in Fig. 3.18 b except for FGGESİU pea 
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water. Problems with the FGGE humidity analysis contribute 
to these large residuals. In the ECMVF analysis scheme, 
humidity observations received a very large veight compared 
vith the model first guess during the period of this storm 
(Paegle, 1983). These humidity observations vere comprised 
of Tradiosonde data and surface observations of sometimes 
questionable quality, vhich led to very noisy analyses. 

At later times, a net moisture sink due to precipi- 
tation occurs above a shallov moisture source due to evapo- 
Eu from the sea (Fig. 3.23 b,c). This precipitation is 
associated vith the diabatic heating in Fig. 3.22 a. As for 
the NOGAPS case, most of the vater vapor accumulation comes 
from horizontal flux convergence near the surface, and from 
1 "£ical flux convergence above about 900 mb (Fig. Jm 3 
277 Eddy terms (Fig. 3.23 h,i) are small. Precipitation 
Busunts are about 30% of those for the NOGAPS case, even 
Bosch the mass circulations are of Similar strength for the 
two cases. This 1s a result of the much lower temperatures 
associated with the FGGE case. It is also possible that the 
NOGAPS model parameterization results in an overestimation 
of precipitation. 

The thickness tendency analysis for this storm is 
presented in Fig. 3.24. The terms are calculated from (3.4) 
as before, except that the tendency term iS computed 
directly from the gridpoint tendencies, and the diabatic 
term is inferred as a reSidual. At the initial time, advec- 
tive warming occurs to the south of the incipient low, with 
25 Oling to the northvest (Fig. bx 15 - mek. thickness 
Meme@ency (Fig. 3.24 j) shows cooling to the north of the 
disturbance which is consistent with the southward movement 
mee jet Streak and associated baroclinic zone (Figs. 3.5 
F Dd 3.6 a,b). f p p) ete indication of self- 
pment at this time. Adiabatic cooling (Fig. 3.24 d) 
baric heating (Fig. 3.2& g) occur to the north of the 


Mow center. 
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25  . 1?) in Fig. 3°16 except for FGGE Type I. 
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At the time of most rapid development (1200 GMT 27 
January, see Fig. 3.12), a net dipolar pattern of varming to 
the east and cooling to the west is evident (Fig. 3.24 Е): 
This dipole is modulated by the thermal advection field 
(Fig. 3.24 b), and is indicative of self-development, as 
discussed earlier for the NOGAPS case. Diabatic heating 
(Pie 3.24 hn) offset by adiabatic cooling (Fig. 3.124 
occurs to the northeast of the low center. 

At maturity, cooling occurs throughout a broad 
above the low center (Fig. 3.24 1). This tends to produce a 
more vertical structure to the system. The thermal advection 
dipole has weakened considerably due to the much reduced 


baroclinity of the system at this time (Fic E 0 7 


E. SUMMARY AND DISCUSSION 


The development of the Type I polar lov has been 
described using thermodynamic and vorticity balance con ami 
erations. The disturbance 1S triggered by low-level baro- 
clinic processes. At genesis, weak warm advection associated 
with the initial baroclinic Surface disturbance Tresu1 7777 
ascent ahead of the low. Associated latent heat release 
further contributes to ascent. Surface convergence leads to 
development of the surface low, which further strengthens 
low-level thermal advection. Thus, initial development 
the polar low is a cooperative interaction between the 
growing vertical mass circulation and the thermody mo meis 
processes. Since surface winds are light at early stap c cmimi 
polar low formation, the tendency of the surface flu ccm 
oppose development is small. 

AS development occurs, thermal advection, adiabar 
cooling and latent heating combine to produce hydro r G 
temperature changes which favor the formation of an anə 


fying upper-level geopotential wave. This wave is coupled 
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vith the developing lov-level cyclone. Thus, the sequence of 
self-development, vhich has been previously applied to 
larger scale mid-latitude extratropical cyclones, appears to 
apply on the smaller polar low scale. As self-development 
proceeds, vorticity advection becomes the leading term in 
the vorticity balance. The required increase in vorticity 
generation is accomplished by expansion of the upper trough, 
and the: movement of the polar low to the cyclonic side of 
the jet. 

During the final stages of development, the nonlinear 
Mi ects of vorticity and thermal advection appear to 
generate conditions favorable for later secondary develop- 
ments. Thermal advection results in progressive amplifica- 
tion of the thermal waves which contributes to frontogenesis 
in an area to the southwest of the low center and an area to 
the east. These areas are the expected locations of cold 
and warm fronts. Furthermore, hydrostatic geopotential 
changes r occurs rimarilv fos a result of advective 
thermal changes result in cyclogenesis aloft near these two 
regions. The PVA ahead of these two lobes of cyclonic 
vorticity 1S associated with upper-level divergence and 
vertical ascent, which will contribute to the future devel- 
opment of baroclinic disturbances near these locations. 

Since the Type I low forms on the anticyclonic side of 
the jet, the upper-level divergent flow is initially ineffi- 
ment in producing vorticity changes aloft. Pings inftfal 
BEN Ch of the mass circulation is accompanied by only modest 
changes in the upper-level vorticity structure. This is 
nur to a cyclone described by Pagnotti and Bosart (1985) 
which formed beneath anticyclonic upper-level flow on the 
equatorward side of the jet streak, and was accompanied by 
25 SDE precipitation. They found that the vertical 
Meeent was forced predominantly by latent heating and 
Emermal advection, with very little initial development of 


upper-level PVA. 
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Because the Type I polar low forms in the absence of PVA 
aloft, it is initially detected as a surface disturbance. In 
the NOGAPS case, a pre-existing lov-level vorticity maximum 
provides the initial disturbance. In the FGGE case, the 
origin of the initial perturbation is unclear. Orograp nikis 
(lee-side) effects associated with Greenland may be impor- 
tant. Strong surface heating as cold air moves off ice- 
` covered terrain would also lead to falling surfZace 
geopotentials. However, understanding the deformation fields 
within the parent cyclone appears to be the key to 
successful prediction of Type I polar low cyclogenesis. 

General synoptic features of the simulated polar low are 
consistent with the FGGE analyses (and earlier observational 
studies of maritime cyclogenesis). This demonstrates the 


realism of the model simulations and justifies the detailed 


analysis of model output. Because the model-generated data 
set is complete to computer accuracy, vorticity, heat; 
moisture budget terms can be diagnosed with confidence. In 


comparison, the FGGE analyses are rather noisy, as evidenced 
by the inability to satisfy the fundamental conserva 
equations. One source of uncertainty is that the archived 
temperature data are from the ECMWF forecast model, which 
was used for the first guess. The humidity observations, 
which were of poor quality, were given too much weight in 
the ECMWF analysis scheme during the period of the FGGE 
storms Im addit ion; the FGGE analyses appear to be espe- 


cially deficient near steep topography such as Greenland. 


in summary, it is evident that the Type 1 polar lov 
grows through ур тә ас processes. Locally reduced 
Richardson number in the vicinity of the disturbance pasa 


tive correlation between vertical motion and temperature, 
and the decreasing baroclinity which accompanies) ey aus 
genesis appear to support this conelusiona In the cases 


presented here, the initial formation of the SUr RERS 
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disturbance is favored by reduced static stability and 
enhanced baroclinity, and occurs without the presence of 
Significant upper-level PVA. The formation of the Type I 
polar low on the equatorward side of the upper-level jet 
appears to contradict the conclusion of Reed (1979) and 
Mallen (1979) that polar lows form on the cyclonic shear 
side of the jet. The complete record from the model simula- 
tion has permitted a detailed description of the Type I 
polar low, nich does not initially have significant PVA 
ASE. À similar description of a Type II polar low which is 


triggered by PVA aloft will be given in the next section. 


Os 


IV. CYCLOGENESIS UNDER DISTURBED RUPEE n) 

As described in the previous chapter, the Type I polar 
low forms under undisturbed upper-level flow on the equator- 
ward side of the upper jet. Development iS initiated by 
low-level baroclinic processes which couple to devela) 
upper waves. Storm intensification results from a coopera- 
tive interaction between the increasing surface circulation 
and the amplifying upper-level features. 

Type II polar lows typically form on the cyclonic ии 
of the upper-level jet in a strong polar air stream, and are 
distinguished from Type I disturbances by the initial pres- 
ence of appreciable upper-level PVA. The polar lows 
described by Mullen (1979, 1983) and Reed (1979) fulfill 
these criteria. These North Pacific disturbances form ın272 
polar airstreams behind or poleward of major frontal bands. 
They are associated with strong baroclinity throughout777- 
troposphere, and originate as areas of locally enhanced 
Convection. A substantial depth of the troposphere Es 
conditionally unstable and strongly heated from below, espe- 
cially at early stages of development. As these storms 
develop, they acquire a comma-shaped cloud pattern, and are 
accompanied by a westward-tilting geopotential perturbation 
throughout the entire depth of the troposphere. 

As with Type I disturbances, the Type Il polar lows 
attain considerable intensity and replace the primary low as 
the main circulation center. Occasionally, a wave devel 
on the polar front just ahead of the comma cloud 171520 
low). If the distance between the comma cloud and the 
frontal cloud band becomes small enough, the two 7467700 
combine to form a cloud image which resembles a mature 


occluded, cyclone This is the so-called instant occ mM EEE 
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process that was first described by Anderson et al., (1969). 
Sometimes, the polar low develops occlusion-like features 
Beparate from the main front (Locatelli et al., 7. 
Hüllen (1983) described examples of Trapidly deepening 
cyclones in polar air streams which formed downstream of a 
mobile upper-level trough within or poleward of the upper- 
level jet. These North Pacific cyclones satisfied Sanders 


1 


and Gyakum"s (1980) criterion for "explosive deepening” 
(average central pressure falls in excess of 1 mb/h for 24 
n). 

Two examples of Type II developments will be presented 
Mae this chapter -one from the NOGAPS simulations by 
Sandgathe (1981) and the other from the FGGE data. Because a 
detailed description of the quasi-Lagrangian analyses was 
presented in the previous chapter, comparatively brief 


results will be presented for these cases. 


SYNOPTIC OVERVIEW 
i) NOGAPS Case 


The model-generated surface polar low formed ina 
27” Initially undisturbed, cyclonically-curved northwes- 
terly air stream behind 4000010 clone. (Exe. tel aj... This 
25107 disturbance grew as an intensifying trough within this 
flow and eventually became a closed circulation and replaced 
the parent low as the main low center (Fig. 4.1 b,c). The 
polar low formed beneath the PVA ahead of a developing 500 
Memerough (Fig. 4.1 d-f). The upper trough attained consid- 
erable amplitude by the end of the period (Fig. 4.1 f) with 
little development of the downstream ridge. Except for the 
ler horizontal scale of this system, this development 
appears very similar to the Petterssen (1956) description of 


Byeroezenesis which occurs in reponse to PVA forcing aloft. 
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This cyclone appears to be more vigorous than the 
NOGAPS Type I low. Whereas the Type I disturbance forms in 
an area of light surface wind, the Type II storm forms 
152 -in a surface flow in excess of 15 m s-". A OD. 
wind speeds in excess of 65 m s-! occur above the low center 
(Fig. 4.2 a). The surface pressure at the disturbance center 
fell from 1010 mb to less than 981 mb over a 36 h period. 
Ziring the final 12 h, an 11 mb fall occurred. This deep- 
ening rate is not far short of Sanders and Gyakum's (1980) 
criterion for explosive cyclogenesis. 

Vertical motion fields (Fig. 4.2 d-f) for the model- 
generated polar lov suggest development of the classical 
comma ascent signature. The similarity betveen these model- 
derived ascent fields and the usual satellite-observed 
Exbratropical cyclone cloud patterns is striking. Realistic 
features such as the appearance of a subsiding air slot and 
cyclonic rotation of the cloud system are reproduced. AS 
for the Type I case, the strongest ascent occurs in the head 
Ehe comma, which is located within 100 km downstream of 
Hu surface vorticity maximum. For this disturbance, the 
Minimum surface pressure tends to be located about 1-2 
degrees of latitude poleward of the maximum surface 
25 -ity (Fig. 4.1 b,c). Therefore, the minimum surface 


pressure occurs on the upstream polevard side of the comma 


head, as evident in Fig. A.2 d-f. Bbis conclusion also 
applies to the NOGAPS Type I polar low (see Fig. 3.2 d-f). 
Mee omma tail appears to lie within the surface trough, аз 


Eo noted by Reed (1979). 
Noe 0Dubesobvpachandatysis (Fig. 4.2 a-c) shows 


En Ucrhe cyclone forms in the exit region of a strong (65 


Meee) Advancing upper jet. As development of the surface 
En occurs, the jet becomes located about 500 km to the 
2” Of the surface cyclone (Fig. 4.2 c). The strongest 
2” -city in Fig. 4.l d-f is found in the zone of enhanced 
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shear vorticity immediately to the polevard side of the vind 
maximum (Fig. 4.2 a-c). Huoner.ase to the Type I polar 
lov, vhich spends a significant fraction of the early devel- 
opment stage on the equatorvard side of the upper jet, the 
Type II disturbance forms and develops directly underneath 
the poleward side of the jet. 

The strength and movement of the upper jet agree 
reasonably with Mullen (1979), who composited 22 cases of 
North Pacific polar lows. At genesis, the composite distur- 
bance was located 500 km to the north of the jet stream, 
which had a mean speed of 65 m s-?. At maturity, the mean 
center was located 750 km to the north of a weaker (55 
m S-!) jet streak. Mullen defined genesis as the stage prior 
to the appearance of a comma cloud configuration, when a 


E area of enhanced convection was evident in the 


infrared (IR) satellite imagery. Maturity occurred when the 
cloud pattern was  comma-shaped, but without excessive 
Spiralling. 


The main disparity between the NOGAPS Type II low 
and the observational cases of Mullen (1979) is the location 
of the upper-level jet at polar low genesis. Tr Cud Pen s 
composite, the jet lies 500 km to the south of the incipient 
disturbance. In the NOGAPS case, Bem Low forms directly 
under the upper-level wind maximum. It is possible that 
genesis in Mullen's cases is detected much later than for 
the NOGAPS low, Since considerable vertical ascent is 
already occurring at genesis in Mullen's cases. At later 
Stages, the model jet coincides with the subsident slot, 
m he tail lying to the south of the jet. This agrees 
femmes with Carr and Millard (1985), who composited the 
large scale environment for 68 comma-cloud systems over the 
Great Plains. They found evidence of two upper-tropospheric 
wind maxima -one aligned along the dry slot axis and the 


Amer cutting across the southern section of the comma tail. 
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Hodur (1984) initialized the Navy Operational 
Regional Atmospheric prediction System (NORAPS) with the 
composited data of Mullen (1979) that was typical of condi- 
tions during the genesis stage of North Pacific polar lows. 
With the full diabatic version of the model, he found that 
the simulated polar low developed about 7/750 km north of the 
jet, in agreement with Mullen's (1979) observations. When 
the adiabatic version of the model vas used, the simulated 
polar low formed 1000 km north of the jet. This suggests 
that latent heating may account for development closer to 
the jet. In Mullen's data, the low formed with a mean 
surface temperature of 7°C while the NOGAPS Type II devel- 
oped with a surface temperature of 14°C. Therefore, the 
effect of latent heating would be greater in the warmer 
NOGAPS storm, which may account for the development closer 
GO Ene Jet. 

Thermal fields for the NOGAPS case (Fig. 4.3) show 
strong initial baroclinity, especially at 500 mb (F1£90 72 
dt Development of the cyclone is accompanied by rapid 
amplification of the thermal wave, which results in a strong 
thermal ridge (trough) ahead of (behind) the surface low 
center. At later stages, the thermal patterns are quite 
Similar to the classical frontal cyclone, with enr || 
baroclinity to the southvest and to the east of the sür 7 


center (Pig G) Denen Ји 
2.  FGGE Case 


The Type II FGGE case chosen for this study fommea 
at 1200 GMT 27 December 1978 in the northwesterly air stream 
of a mature 972 mb cyclone which vas located over sa 
central North Pacific Ocean. Synoptic features of Chro cM 
lov are similar to those for the NOGAPS case prevrouza 
described. The surface disturbance developed as a trough in 


the polar air stream (Fig. 4.4 a) and migrated about wea 
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parent low to eventually become a closed circulation (Fig. 
4.4 b,c) in the southwestern Gulf of Alaska. Associated 
surface pressure falls were small, but the relative 
vorticity increased from around 4 x 10-* s=* to 8 700072 
at the low center. The surface trough in the earliest anal- 
ysis (Fig. 4.4 a) is stronger than the corresponding initial 
NOGAPS disturbance (Fig. 4.1 a). In the NOGAPS case, the 
3-hourly data permitted detection of the incipient polar low 
at a much earlier stage than for the 12-hourly FGGE data. 

As inthe NOGAPS case, the surface disturbance 
formed beneath the area of PVA ahead of a rapidly developing 
500 mb trough (Fig. 4.4 d). At the end of the period, an 
almost-closed upper-level low was present directly above the 
surface lov (Fig. 4.4 f). This vertically stacked orris 
tion of the lower and upper features was not evident for the 
NOGAPS storm. In that case, the 500 mb low center 778 
located 800 km to the northwest of the surface low 45502 
final stages of development (Fig. 4.1 f). 

Surface and 500 mb thermal fields (Fig. 4.5) showed 
that the polar low formed within a deep baroclinrc zone mM 
convoluted patterns of surface temperature at IS8ter jsp || 


of development (Fig. 4.5 b,c) are probably due to data 


uncertainty. Since the model-zenerated surface temperature 
fields in Fig. 4.3 a-c are known to computer accuracy msn a 
are smooth in comparison to the observed temperatures. At 


the initial stage, a cold pool with temperatures below -32 6 
was evident at 500 mb. Based on Fig. 4.5 a,d, the St EEs 
stability index at this time vas 2.4”K per 1005507022 The 
moist adiabatic stability index corresponding to the) Sun ae 
temperature of 9°C was 3.5"K per 100R TF which indicated 
that the incipient polar low environment was conditiSo wm 
unstable. Later, the 500 mb cold center was located 52207 


south of the surface low center (fie. 777 
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The isotach analysis (Fig. 4.6 a-c) indicates that 
the polar lov formed at the left exit region of a 50 m s-" 
2” streak at 300 mb (Fig. 4.6 a). After 24 h, a 65 m s-? 
jet streak vas located about 1000 km to the south of the 
BHusturbance center (Fig. 4.6 b). The strength and position 
of the upper-level wind speed maximum is similar to that for 
the NOGAPS Type II case. The maximum wind speed aloft is 
nearly double that for the two Type I polar lovs studied. 

The ascent which occurred to the east of the surface 
development (Fig. 4.6 d-f) does not exhibit the comma-shaped 
signature shovn by the NOGAPS case (Fig. 4.2 d-f). The more 
2  :P€6E analysis is probably due to data uncertainty. 


Upward motion is generally weaker than for the NOGAPS storm. 


Summary 


The Type II polar low forms in strong surface flow 
on the cyclonic side of the upper jet, and is accompanied by 
appreciable baroclinity throughout a deep layer. The use of 
model-generated data permits early detection of the inci- 
pient surface wave and the accompanying upper-level flow, 
and shows that a pre-existing upper-level trough and associ- 
ated PVA are required for formation of the Type II distur- 
bance in the model. Unlike the Type 1 lovs described 
earlier, there is no initial closed circulation associated 
with the Type II surface disturbance. Thus, it may be diffi- 
cult to detect the early stages of such a development with 
observational data. -At later stages, it is not possible to 
discriminate between Types I and II, since an amplifying 
upper-level wave and associated PVA is coupled with the 
intensifying low-level cyclone in both cases. 

A major difference between these Type II distur- 
Dances and the Type I polar lows studied earlier is the 
Stronger initial flow in which the Type II low is embedded. 


Wind speeds vary from around 15-20 m s-! at the surface to 
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50-60 m s-! at 300 mb. Corresponding typical values for the 
Type I polar lows are 0-10 m s-' at the surface to around 35 
mes-' aloft. Although vertical wind shear (and therefore 
baroclinity) is similar for the two polar low types, the 
initial presence of strong surface wind suggests that 
surface fluxes of heat, moisture and momentum (friction) are 
important in the formation of the Type II polar low. 

In the next sections, mass, vorticity, heat and 
moisture budgets will be presented to illustrate development 
of the Type II polar lovs. An additional obiective of this 
analysis is to determine the reasons vhy some polar lovs 
form under straight upper flov (Type 1) vhile others require 


the presence of upper-level PVA at genesis (Type 11). 


PBB MASS AND VORTICITY ANALYSIS 
1. NOGAPS Case 


Mass transport terms for the simulated polar lov are 
presented in Fig. 4.7. The horizontal mass flux convergence 
mame 74.7 a) has a two-layer structure similar to the Type I 
disturbance (Fig. 3.7 a) with a deepening convergent layer 
below a layer of mass outflow. The kinematically-derived 
Mertical velocities (Fig. 4.7 b) increase strongly during 
the first nine hours of development, and are accompanied by 
increasing low-level convergence (Fig. 4.7 a). This is 
followed by a period of weaker, but pulsating growth. Thus, 
Eeowtlh of the mass circulation is more rapid than for the 
Type I case. As described earlier, the 12-h oscillation in 
20 ass circulation is an artifact of the external gravity 
wave propagating within the NOGAPS model. 

Cyclogenesis commences near the surface and later 
ENuds throughout the troposphere (Fig. 4.8). Initially, 
there is weak cuu luconesis alotfE? Vorticity increases in 


the lov troposphere coincide vith enhanced low-level 
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convergence (Fig. 4.7 a). Thus, development of the surface 
cyclone appears to occur in response to the increasing mass 


circulation evident in Fig. 4.7. 
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Fig. 4.8 vorticity tendency (10-*% s-!) for NOGAPS 


Type II 


The peaks of low-level cyclogenesis at 168 h and at 
mae (Fig. 4.8) correlate well with the mass circulation 
Maxima (Fig. 4.7) at these times. However, this “pulsing” of 
vorticity increases appears to be partially modulated by the 
meee epravity wave which is an artifact of NOGAPS. «S 
illustrates a disadvantage of the 3-h data resolution. At 
mes resolution, much of the inherent noisiness of the model 
becomes evident. However, time averaging can be used to 
52566r out this undesirable noise, vhile preserving the 
slover meteorological modes. 

Vorticity budget results are presented in Fig. 4.9. 
Over the first 12 h, the upper-level PVA term (Fig. 4.9 d) 


exceeds 10 x 10-? s-! per day. It is this initial presence 
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of strong PVA which characterizes the Type II polar low. By 
comparison, the average PVA for the Type I storm is less 
than 4 x 10-? s-! per day over the first 21 h of development 
(Fig. 3.10 d). PVA for the Type I storm becomes large Joma 
at later stages of development as amplification of the upper 
waves occurs in conjunction with self-development. 

The large PVA at early stages of storm development 
is offset by strong upper-level divergence (Fig. 4.9 aim 
This is consistent with the growth of the vertical ies 
Circulation düringe mes and the resulting generation 
of low-level vorticity. Twisting (Fig. 4.9 g) accounts ai omaa 
significant loss of vorticity below 500 mb. This 15 ares 
of the localized nature of the initial vertical morron 
(large horizontal gradients of omega) ina strongly baro- 
Clinic environment? 

As development occurs, vorticity increases oc m 
throughout the depth of the troposphere (Fig. 4.9 07702 
Aloft, this generation is largely the result of PVA, and is 
partially offset by divergence (Fig. 4.9 e,f). Nearer the 
Surface, generation is due to convergence with sinks arising 
from vertical eddy divergence (Fig. 4.9 h,i) and fric mus 
Thus, the vorticity budget at later stages of this case is 
very similar to that for the Type I storm (Fir = FIER 

The horizontal distributions of the terms ın 2 
vorticity equation at 500 mb are shown in Fig. 4.100770772 
initial time, PVA (Fig. 4.10 a) dominates the area direz 
above and to the vest of the surface disturbance. This 45 
compatible with a strong vorticity gradient alone then riks 
in these locations (Fig. 4.1 d). Since the wind veloc rromi 
this level exceeds the translation speed of the disturbance, 
the quasi-Lagrangian vorticity advection is Targe and?2722 
tive. The maximum contribution from the divergence term 
(Pigs 4A 10000 occurs a little to the east of the cere = 


Vertical advection and twisting effects (not shown) are 
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small. The net quasi-Lagrangian vorticity tendency (Fig. 
4.10 g) is an increasing (decreasing) vorticity to the west 
(east) of the center, vhich indicates amplification 67207 
,upper wave. This dipole of vorticity tendency 1S conside 
ably larger than that at the initial stage of the Type I 
polar low (Fig. 3.11 g). This suggests that the amplia 
tion of the upper wave is more rapid for the Type 11 polar 
lov. 

At later stages in this storm, strong PVA extends in 
a band from directly above the surface low center to east of 
the low (Fig. 4.10 b,c). This is consistent with the east- 
ward extension of the area of positive vorticity in Fig. 4.1 
e. The strong PVA to the southwest of the lov center is 
associated with a mobile short wave trough just oft fM 
western edge of the domain. The vigorous mass circulation 
causes a strong dipolar contribution from the divergence 
term (Fig. 4.10 e). In the trough, convergence is very effi- 
cient in generating vorticity due to the already high values 
of vorticity (see Fig. 4.1 e). Although the divergence sink 
exceeds PVA ahead of the storm, an additional posu 
contribution from both vertical advection and tvistine 6 
shown) occurs in this location, and results in a nera 
tive tendency (Fig. 4.10 h). The tendency patterns at me MES 
time, and at the end of the sequence (Fig. 4.10 i) are 
remarkably similar to those at the end of the Type 1 NOGAPS 
storm (Fig. 3.11 i), with positive tendency to the ed oE M 
southwest of the center. This suggests that the later evolu- 


tions of the Type I and II storms are s mi 
2. YFGGE Case 


The mass budget analysis for the FGGE case shovs a 
deepening layer of convergence below upper-level outflow 
(Fig. 4.11 a). The strength of the mass circulation pela 
1200 GMT 28 December, and decreases thereafter (Fig. 4.11 
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b). The minimum area-averaged omega of just under -5 mb per 
h is about half that for the NOGAPS case. In the NOGAPS 
case, a much broader maximum in the mass circulation 
evident. Omega values below -8 mb per hour extend over a 
30-h period (Fig. 4.7 b). A similar difference was observed 
between the FGGE and the NOGAPS Type I cases described 
earlier (compare Figs. 3.7 and 076 This difference high- 
lights a model weakness. In nature, strong vertical ascent 
associated with a particular weather system is usuaily 
Short-lived as the system undergoes the cycle of growth, 
maturity and decay. This is evidently the case for the FGGE 
examples (see Figs. 3.8 and 4.11). In contrast, the NOGAPS 
model appears to systematically sustain strong vertical 
ascent over a period of 2 to 3 days (Figs. 3.7 and aR 
Another possible explanation for the weaker vertical motion 
in the FGGE case is that the vertical motion is averaged 
over a larger area due to data distribution. 

The vorticity tendency time section for this Type 1l 
FGGE case (Fig. 4.12) s ows that cyclogenesis commences near 
the surface and extends throughout most of the troposphere 
by the end of the period. This cyclogenesis peaks at 0000 
GMT 29 December, which is about 12 h after the maximum 
the mass circulation. The strongest vorticity tendency is 
about half that of the more vigorous NOGAPS Type 119260 
(see Fig. 4.8). 


Vorticity budget results are shown in Ејр MAMMA 
During the first 12 h of development, a shallov surface- 
based layer of cyclogenesis occurs (Fig. 4.13 a, see main 
Fig. 4.12) as a result of low-level convergence (FIg NE 
d). This vorticity generation is strongly offset biuk a 
tion, as inferred from the negative vorticity budget resi- 
dualumu TES 4.13 a. In most bulk representations r 
PRES friction appears in the momentum equation as the 
vertical gradient of the vector eddy stress, which in 07 
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is proportional to the product of the Surface wind Speed wag 
the surface wind vector. As a first approximation, the iS 
tional sink of “vorticity is proportion EE product of 
the surface vorticity and the surface wind speed. Thu?” 
frictional dissipation of cyclonic vorticity is larges 2 
the surface flow is strong and cyclonic, as is the cas D 
Fig. 4.4 a. Friction opposes development of a disturbance in 
the strong cyclonic surface flow such as in the Pi r 
stage of Type II polar lows. By contrast, surface ftri m 
is minimal at initial stages of the Type I polar low, 3007 
surface wind speeds are small. 

The vorticity budget at early stages of the NO TIHE 


Type II polar low indicates a small positive residual (Pages 


4.9 a), which would appear to contradict this arzgzumanaş 
However, in the NOGAPS case, the initial area-averaged 
vorticity is lose e Te o It is possible that frict Tom 
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dissipation of cyclonic vorticity near the 3727707772 
center is more than offset by dissipation of anticyclonue 
vorticity in surrounding areas included in the budget 
volume, which results ina small positive budget residual. 
In that case, and in both the Type I vorticity budeeme 
(Figs. 3.10 a-c and 3.21 a-c), the budget residual becomes 
increasingly negative at later stages of development. This 
indicates increasing frictional dissipation Of Cyc Omis 
vorticity as the surface circulation incre Ice of 

Aloft, the FGGE vorticity budget indicates ada 21722 
vorticity sink from divergence, which is offset by PV SR DE 
4.13 d,e). At later stages, the contribution from the 
gence term is smaller due to the reduced vertical mass 
circulation (see Fig. 4.11). The vorticity budgets for 75 


case, and for the Type I FGGE case (Fig. 3.21) yield resi- 


quals com comparable magnitude to the leading terms. 
Therefore, the estimates of the various budget terms are of 
ФипБтопбитассатасы By contrast, the model-generated data 


permitted reliable diagnoses of vorticity budget terms. 


7 isetssi6u 
The evolution of- vorticity EDIT has been 
described for two Type II polar lows. The analysis shows 


that PVA associated with an advancing upper-level trough is 
important at early stages of development. Апр1:Ејсастевини 
the upper wave is evident at early stages, especially inm the 
trough where convergence 1S an -efficient source of 
vorticity: Im contrast, the Type I polar low LormS UNAS 
straight or anticyclonic upper flow with initial ps 5 
the downstream ridge. Because the vorticity in the location 
of the upper trough iS comparatively low for the IPM 
polar low (see Fig. 3.1d), convergence is less effective in 
amplifying the initial trough. In the Type 1 case, the ua 


ridge andthe divergence are co-located, so the tae 
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intensifies rapidly at early stages, in spite of the lov 
vorticity. In the Type II case, the initial presence of PVA 
offsets the divergence sink, and development of the dovn- 
stream ridge is minimal. At later stages, evolution of the 
upper waves progresses similarly for both types. Finally, 
surface friction appears to strongly oppose development of 
the Type II disturbance due to the strong surface flow. The 
presence of this strong surface flow represents an important 


difference between Type I and II polar lovs. 


C. THERMODYNAMIC ANALYSIS 
1. NOGAPS Case 


Results of the heat budget calculation are shown in 
Fig. 4.14. At initial stages of development, potential 
temperature rises occur throughout a deep layer (Fig. 4.14 
a). Mos wimi is hydrostatieally consistent with the 
weal vorticity tendency profile (Fig. 4.9 a), which has 
the strongest  cyclogenesis near the surface. The largest 
EuwErrbution to this warming is from the diabatic heating, 
which is offset by the mean and eddy flux divergence terms 
Mere. 4.14 d). The diabatic heating at initial stages of 
5:0 ipe I polar low (Fig. 3.15 d) is small by comparison. 
This reflects the more rapid growth of the vertical mass 
circulation in the Type II case. The budget results at 
later stages of development are very similar to those for 
the Type I low (see Fig. 7777777 176 leading terms at this 
stage are diabatic heating and the offsetting mean heat flux 
convergence (Fig. 4.14 e,f). The eddy term (Fig. 4.14 h,i) 
is small because the budget volume includes both cold and 
Warm advection. 

The heat budget shows that the most significant 
differences between Type I and II polar lows occur at early 


stages of development. The diabatic heating for the Type II 
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Polar low is significant at genesis and reaches a maximum 
during the first 6-8 h of development (Fig. 4.15 a). By 
contrast, much slover grovth of diabatic heating is exhib- 
ited by the Type I storm (Fig. 3.18 a). At later stages of 
development, evolution of heat budget terms is similar for 
the two polar low types. 

The Type II NOGAPS case appears to develop more 
rapidly at early stages than the Type I case. The early 
peaks in cyclogenesis (Fig. 4.8), diabatic heating (Fig. 
a) Sand mass circulation (Fig. 4.7 b) correlate well 
with the strongest baroclinic conversion (Fig. 4.15 b). A 
Mater maximum of baroclinic conversion coincides with the 
second "pulse" of cyclogenesis. This correlation between 
storm intensification and the baroclinic conversion suggests 
that the storm grows as a baroclinically unstable wave. 

A notable feature of the diabatic heating distribu- 
won in Fig. 4.15'a is the large surface heating rate of 7?K 
per day at polar low genesis. Large surface heat fluxes 
(Fig. 4.16 a) are caused by high surface winds (in excess of 
Nm s-* near disturbance) moving over warmer water (see 
mie. 4.3 a). Since the fluxes are largest in the coldest 
air, their effect is to reduce the horizontal temperature 
— nt in Fig. 4.3 a. At later times, a dipolar pattern of 
surface heating  (cooling)(Fig. 4.16 b,c) develops behind 
(ahead of) the intensifying surface trough. This heating 
pattern is in antiphase with the amplifying surface tempera- 
ture wave (Fig. 4.3 a-c) and therefore opposes development. 
By comparison, surface fluxes for the Type I polar lows are 
weak because of light surface winds, especially at early 
Stages of development. The presence of large surface 
2 nç is an important characteristic of the Type II polar 
low during formation and early development. 

The moisture budget (Fig. 4.17) 2555 6565 that large 


Meeumulations of water vapor occur primarily as a result of 
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the convergence of horizontal and vertical mean a. (Fig. 
EU d-f). The eddy terms (Fig. 4.17 g-i) are’ weak by 
comparison. As this accumulation of water vapor exceeds the 
increase in vapor content, cM Coup Sink (precipitation) 
exists in low-mid layers (Fig. 20757: At the initial 
stage, the mean vapor content of the budget volume is 
increasing (Fig. 4.1/ a). This is consistent with the 
initial warming trend (Fig. 4.14 a), since the Saturation 
specific humidity increases with temperature. Therefore, a 
Smet icant fraction of the flux convergence of water vapor 
Semeributes to an increase in the moisture content. The 
shallov surface layer of positive residual (Fig. 4.17 a) 
indicates that evaporation from the ocean is responsible for 
some of this moistening. Above the surface layer, the 
effects of horizontal convergence at lower levels, and 


E rcal flux convergence at higher levels (Fig. 4.17 d) are 


L9 


sufficient to ensure continued pPprecipi2:57 77 These flux 
convergences in turn are the result of the increasing 
vertical circulation. The surface evaporation ensures that 
the low-level inflow is sufficiently moist to sus 7 
continued precipitation: 

The terms in the thickness tendency equation (3.4) 
have been computed for the model-generated Type II distur- 
bance (Fig. auus Quasi-Lagrangian thermal  advection 
patterns at the initial time (Fig. 4.18 a) present a dipole 
of warming ahead of and cooling behind the developing 
cyclone. Although this seems to contradict Fig. 4.3 2390 D 
indicates cold advection in all sectors), the thermal advec- 
tion in Fig. 4.18 ais calculated relative to the SON 
which is translating to the east at about 23 m s- Wk n 
this translation velocity is subtracted from the winds, (ene 
storm relative motion ahead of (behind) the cyclone 7 
toward the cold (varm) ut. The intense warm advection to 
the northwest of the center is also a cənsequence of the 
qudası Lagrangian calculdueron The diabat c heating to the 
east of the lov center (Fig. 4.18 g) almost exactly cancels 
the adiabatic cooling (Fig. 4.18 d). As a consequence, the 
pattern of net thickness changes (Ereg ^8 gm closely 
matches the thermal advection (Fig. 4.18 a). Since cool 


(warming) of a column results in falling (rising) gcor EE 


tial at they top: Gf thes comunn- these patterns result in 
amplification of the upper waves. This amplification mM 
evident in Fig. 4.10 g as discussed earlier. Thus, 57 


contributions of thermal advection and diabatic effects 
self-development as described earlier for the Type I polar 
low also apply to the Type IL diSsturbaner” 

These self-development tendencies occur in a favor- 
able environment of low static stability. Averaged over the 
budget volume, the stability index is 3.7?*K per 1002552 


genesis, which is very close to the moist adiabatic aoe 
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Adiabatic warming associated with subsiding air near the 
vestern boundary of Fig. 4.18 d opposes the cold advection. 
Because this area is in the cold air vhere surface fluxes 
are large (see Fig. 4.16 a), the static stability is reduced 
and the braking effect due to subsidence is also reduced. 

At later stages of development, there is a cyclonic 
Eotatron of the thermal advection patterns. This rotation is 
associated with amplification of the thermal wave. The net 
thickness tendency pattern (Fig. 4.18 k) is strongly fronto- 
lytic in the vicinity of the storm center, and frontogenetic 
to the south and southwest. This is consistent with the 
southward displacement of the upper jet streak relative to 
the center, Since the largest thermal and geopotential 
gradients occur with the jet. Since the 500 mb flow above 
the low is southwesterly at this time, the thickness changes 
in Fig. 4.18 Kk indicate continued self-development of the 
upper-level features. 

Anal stages of the storm, the net tendency 
dipole (Fig. 4.18 1) has rotated further with largest 
cooling to the southeast of the center and weak warming to 
the northwest. This tends to produce a more vertical geopo- 
tential structure. Averaged over the domain, net cooling 
Eis. in accordance with the heat budget (Fig. 4.14 c). 

Cross-sections of potential temperature and 
Richardson number through the storm center at genesis are 


presented in Fig. 4.19. The reduced static stability in the 


western sector of Fig. 4.19 a, as implied by . the widely 
Spaced potential temperature contours, is consistent with 
the area of low Ri to the west of the low center. This is 


consistent with the preferred development of the upper 


25060 and is a result of the strong surface heating in this 


EEGror. The deep layer of low Ri to the north of the low 
Eon (Fig. 4.19 d) is a combination of low static 
Seaotlity below about 600 mb, “nü ennanced baroclinity 


J 


aloft. The jet is located directly above the surface iiis 
position at this time. Therefore, the favored location e 
baroclinic development is in the reduced Ri. layer just pole- 


ward of the jet. 
2. FGGE Case 


Heat budget results for the  FGGE case are presented 
in Fig. 070 As with the Type I FGGE heat budget, the 
residual (Fig. 4.20 a-c) includes the effects of din m 
heating as well as data uncertainty and interpolation error. 
The horizontal eddy term (Fig. 4.20 g-i) is small because 
the translating budget area includes both cold and warm 
advection. The large contribution in the upper troposphere 
from the net mean flux convergence (Fig. 4.20 e,f) is not 
balanced by any other budget term. Since the mass budget is 
balanced, uncertainty in the temperature data may contribute 
to the anomalous magnitude of this term. In the NOGAPS heat 
budgets, a similar large negative contribution from tie mean 
fluxes is offset by strong warm advection aloft nE 5 
7506775) 

The time evolution of the heat budget residual term 
1S presented in Fig. 4.21 a. Occurrence of diabatic heating 


in the upper troposphere of the intensity implied in Fig. 


& 21 a as unrealistic. Therefore, this residual reflects 
data and computational uncertainty. In the heat биа. 
involving the complete model-generated data, considerable 


computational error occurred in the upper troposphere (Fags 
3.15 b,c and 4.14 b,c). This error occurs during 65-77 
polation from sigma to pressure levels, as suggested 
earlier. In the production of the FGGE level III-B lar M 
used “ın a Eh S es ue n the first-guess predictions were 
obtained in sigma coordinates from the ECMWF numerical model 
(Bengtsson et al., 1982). Interpolation of these field mEn 


pressure levels used cubic splines, as for the NOGA 
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melds A possible additional source of uncertainty is the 
use of ECMWF model initialized temperatures in this study. 
The departure of these temperatures from layer temperatures 
computed hydrostatically from the geopotential data is 
largest in the stratosphere (Wash and Calland, 1985). This 
temperature uncertainty aloft undoubtedly contributes to the 
large residuals in Fig. 4.21 a. 

At lower levels in the troposphere, diabatic heating 
from latent heat release and surface heat fluxes contributes 
to the heat budget residual in Fig. 4.21 а. In the NOGAPS 


cases, the diabatic heating was strongly modulated by the 
vertical velocity (note the similarity between Figs. 3.7 b 
and 3.18 a, and between Figs. 4.7 b and 4.15 a). The 


dissimilarity between the time sections of the heat budget 
residual (Fig. 4.21 a) and omega (Fig. 4.11 b) for this FGGE 
case shows that the residual is contaminated by data uncer- 
tainty. Strong surface fluxes at early stages possibly 
contribute to the positive residuals near the surface. At 
later stages, the negative residual at low levels appears to 
indicate surface cooling as the storm moves rapidly poleward 


(see Fig. 4.5 c). 


Meba oCclinic conversion term (Fig. 4.21 b) indi- 
7 s a thermally direct circulation (warm air rising, cold 
air sinking) throughout most of the troposphere. This is 


consistent with the previous cases studied and shows that 
mise low grows through baroclinic processes. The magnitude 
of this term is about half that of the more vigorous Type II 
157 PS storm (Fig. 4.15 b). 

The terms in the thickness tendency equation (Fig. 
E indicate the possibility of self-development. Although 
the braking term (Fig. 4.22 d-İf) and the diabatic heating 
200 4-22 g-i) are rather noisy, patterns of thermal advec- 
Es ic 4.22 a-c) and net thickness tendency (Fig. 4.22 


j-1) are more coherent. Tü us3l”feetion contributes 
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warming (cooling) ahead (to the rear) of the surface low, 
which is favorable for self-development. As in previous 
cases examined, these patterns rotate cyclonically with 
time. 

This analysis is subject to considerable uncer- 
tainty. Of the four terms in the thickness tendency equa- 
tion, probably only the thermal advection is derived with 
any degree of confidence. The braking term is the product of 
tvo terms (vertical velocity and static stability) which 
involve horizontal and vertical differentiation and vertical 
integration. The diabatic term is obtained as a residual. 
Bomparison of Figs. 4.22 d and g indicates that diabatic 


heating exceeds adiabatic cooling in the ascent region near 


the low center. Miss cons stent with the conditionally 
unstable lapse rates noted earlier. Only the net thickness 
tendency in Fig. 4.22 k shows thickness changes which are 


indicative of self-development. Because thickness tendencies 
are computed directly, tendencies at the beginning and end 
of the sequence (Fig. 4.22 j,1) are less reliable than those 
during the sequence (Fig. 4.22 k). In both the Type I and 
Type II NOGAPS examples, net thickness changes occurred 
which were in phase with the advective term. Using the 
model results as a guide to the interpretation of the FGGE 
case, it is probable that the sequence of self-development 
also occurred in this observed Type II polar low, as 
suggested by Fig. 4.22 k. 


D. DISCUSSION 


Mire quasi-Lagrangian analysis has shown that the Type II 
polar low is a baroclinic phenomenon initiated by upper- 
level PVA and intensified by low-level thermal processes. 
Continued storm intensification results from 505 


development, as for the Type 1 case. Because of the 


J 
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Fig. 4.22. As in Fig. 4.18 except for ECCE”777 7 
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Mera lily strong surface vinds, surface friction, heat and 
moisture fluxes are important at early stages Ot 
development. 

In spite of the larger translation speeds of the Type II 
storms; the quasi-Lagrangian analyses yield vorticity and 
thickness tendency patterns that are qualitatively similar 
to the slower Type I polar low. In an Eulerian frame, local 
thermal and vorticity changes result from both wave amplifi- 
ENESOn and translation. However, when viewed from a 
quasi-Lagrangian perspective, the large translation contri- 
bution is removed, and wave amplification is similar for the 
two types. Because baroclinic development is dependent on 
vertical wind shear, there is little apparent difference 
between a rapidly moving storm growing in a strong surface 
flow beneath a stronger upper jet (Type II) and a slow- 
moving disturbance in weak surface flow beneath a moderate 
upper jet (Type 1). 

In the maritime atmosphere, Surface fluxes of momentum, 
heat and moisture strongly modify the large-scale environ- 


ment. Because these fluxes are dependent on the surface wind 


velocity, they represent an important difference between 
Type I and II disturbances. Hü eronz cold air outbreaks 
2550 ar to Figs. 4.1 a and 4.4 a, enhanced surface friction 


opposes polar lov development. Strong surface heating on the 
cold side of an intense baroclinic zone tends to reduce the 
surface temperature gradient to a magnitude comparable to 
2 51 gradient. In addition, the heat fluxes oppose the 
thermal advection dipole associated vith a perturbation in 
i Surface flow. Therefore, the surface heat flux tends to 
E23 stabilizing influence on weak perturbations in the 
520 ce flov. This is consistent with Sandgathe (1981), who 
momma that addition of the NOGAPS sensible heat flux weak- 
ened the low-level temperature gradient and resulted in less 


илбиз secondary lows. The damping effects of friction and 


LL 


Surface heating may explain why polar low developments do 
not occur in all strong cold outbreaks. However, Mansfield's 
(1974) conclusion that polar lows develop only when the 
surface winds are light is not supported by the cases 
Studied here. 

The tendency of friction and surface heat fluxes to 
Stabilize the low-level flow suggests why the short-wave 
trough and associated PVA are required for the Type II polar 
low. When a perturbation in a strong surface flow develops 
in spite of the stabilizing influence of the large heat 
fluxes, an upper-level divergence maximum is required to 
sustain a vertical mass circulation: Surface convergence 
induced by this ascent causes low-level vorticity generation 
that leads to formation of a surface trough. às this Tow: 
level circulation increases in strength and vertical extent, 
the growing thermal advection eventually results in the 
onset of self-development. 

A second contribution of the strong surface near 
Suggested when the development occurs in a lov static 
stability environment. An advancing trough aloft (with 
ahead of it), coupled with strong surface heating; Tesis 
ina further reduction of static stability in the Comm ии 
below the trough. Consequently, the subsidence heating which 
would have offset the strong horizontal advection alot mame 
reduced. With larger net cooling aloft, rapid growth Ope 
upper trough occurs. Thus, surface sensible heating may play 
a dual role when the environmental static stabi iC S 
already rather low, as is usually the case with the podem 
low (Mullen, 1979; Reed, 1979). Addition of heat neam 
Surface in antiphase with the growing thermal perturbation 
opposes self-development at low levels. However, strong 
surface heating in a near-neutral environment may lead to 
further reductions in static stability throug mot e 


layer of thé troposphere. For the proper combination gm 
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5 -"ronmental static stability and strong surface heating, 
the self-development of the upper-level trough may be 
enhanced. In these Special circumstances, it 1s hypothesized 
that strong surface heating in the cold air may indirectly 
contribute to surface low development. 

Since formation of the Type II polar low is triggered by 
the vertical mass circulation, the reduced static stability 


region on the cold side of the jet is a favored location for 


development, especially when the surface heating is large. 
For a balanced upper-level jet, divergence is co-located 
with the PVA maximum near the left exit. For these reasons, 


the Type II low forms below an area of PVA on the poleward 
Side of the upper jet. Intensification of this low-level 
cyclone occurs as it becomes coupled to a developing upper- 


level wave. 


In the Type I polar low, the initial surface winds are 
light and the air-sea fluxes are small. In a baroclinic 
environment, the thermal advective effects are not signifi- 


emeliy modulated by the stabilizing effect of the surface 


fluxes and the perturbation grows. Therefore, the distur- 
bance forms near the maximum surface baroclinity. Because 
the baroclinic zone slopes toward the cold air, the surface 


low will form on the equatorward side of the upper jet. 
Since the Type I low forms as a low-level baroclinic distur- 
bance, no initial PVA aloft is required. Subsequent develop- 
ment occurs in conjunction with an amplifying upper-level 
wave which is coupled to the developing surface cyclone. 
Evaporation from the surface is an essential contribu- 
Zu o the moisture budget of the developing polar low. 
These fluxes ensure an adequate supply of moisture for 
pseerpitation. At later stages of development, latent heat 
release is essential for continued development. When 
Sem@eegathe (1981) withheld surface moisture fluxes, the 


secondary lows formed, but failed to develop significantly. 


mn. 


The heat and moisture budgets in this study verify that 
surface moisture flux is essential for polin lov 


development. 
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V. COMPARISONS WITH NON-DEVELOPING CASES 


In Chapters III and IV, it was shown that Type I and II 
polar lows grow through baroclinic processes and are modi- 
fied by surface fluxes and latent heat release. Formation 
and growth of these disturbances is favored by the presence 
of reduced static stability, senonossbaroclinity and. high 
humidity. However, before it can be concluded that these 
conditions are necessary and sufficient for cyclogenesis, 
the hypothesis that the absence of one or more of these 
conditions inhibits  cyclogenesis is examined. Foreprhis 
purpose, two non-developing polar lovs from the model- 
generated data set are examined. 

In nature, there are many non-developing polar lovs that 
remain as secondary disturbances and exhibit pressure defi- 
cits of only a fev mb (Forbes and Lottes, 1982). Such non- 
developing disturbances may initially resemble Type I or II 
polar lows, but do not develop to the same degree. These 
disturbances continue to migrate about the parent low and 
are eventually absorbed into the mean flow. Comparisons 
BX. the initial state of these weakly developing distur- 
bances in the numerical simulations and the more 775255 
developing Type I and II cyclones will help to confirm the 
25:01 various physical processes in the growth of polar 


lows. 


57” 0) DEVELOPING TYPE I DISTURBANCE 


1. synoptic Description 


Selected synoptic fields for a non-developing Type I 
more Ow are shown in Figs. 5.l and 5.2. This disturbance 


Мил 2 Strongly baroclinic surface environment (Fig. 5.1 


Лә 


a) under the right exit region of a 50m sS SEM UN 
(Fig. 5.2 a). The disturbance subsequently becomes embedded 
in the northwesterly flow (Fig. 5.1 bo around the parent 
cyclone, which is situated far to the northeast, and 45 
eventually absorbed into the mean flov (Fig. 5.1 c) vithout 
significant development. The surface disturbance remains 
well to the equatorward side of the 300 mb jet (Fig. 5.2 
b,c) thrcuehourubis £01967 

The initial 500 mb flov (Fig. 5.1 d) is only "eam 
disturbed, vith the highest vorticity to the northvest in a 
broad tongue on the polevard side of the jet streak. NVA 
prevails directly above the lov center, vith an area of veak 
PVA centered about 200 km to the south. This PVA is associ- 
ated vith decreasing shear vorticity on the equatorvard side 
of a weak secondary jet maximum just southeast of the 
disturbance center (rip. xo cM An area of weak upward 
motion (Fig. 5.2 d) occurs in this area. This area of =m 
persists east of the surface center without intensification 
(Fig. 5.2 e,f). Subsequent development of upper features is 
slight, with a weak trough to the west of the disturbamee 
center by the end of the period (Fig. 5.1 f). The areas of 
upper-level cyclonic vorticity and upward motion). near 
eastern boundary of Figs. 5.ld and 2: are associated 
with another more vigorous disturbance wh:ch is just visible 
on the eastern boundary of Fie. 5.1) ae 

Initial conditions for this case are similar In SNE 
respects to the FGGE Type I storm described in Chapter III. 
Strong baroclinity and a location near the right exit Toen 
of an advancing upper-level jet are evident at the genesis 
Stage of the vigorously developing FGGE case (Fig. 3.5 a) 
and the non-developing case described here. However, the 
subsequent evolution of the disturbance in Fig. mS OPES 
markedly from both the FGGE and the NOGAPS Type I cases in 
Chapter III. It is the objective of thus section oK MN 


mine why this disturbance failed to inven ees 
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PINO EE re rey and Heat Budget Analysis 


107767 100 /S1£İüC into factors which hinder develop- 
ment of the Type 1 disturbance described in this section, 
the results of the quasi-Lagrangian budget analysis of this 
case vill be compared vith the more vigorous Type 1 storms 
analyzed in Chapter III. Since reliable estimates of budget 
terms are obtained from the model-generated fields, discrim- 
ination between developing and . non-developing disturbances 
should be possible. Since the FGGE data budgets presented 
earlier were not closed to an acceptably small uncertainty, 
the analysis of non-developing cases vill be confined to the 
NOGAPS simulations. 


The vorticity budget analysis for EhTs non- 
developing Type 1 disturbance shovs a deep layer of 
@ee@measing vorticity (Fig. 5.3 a). Near 300 mb, this is due 


to strong NVA, while a weaker divergence sink 1s evident at 
mid-levels (Fig. 5.3 d). Near the surface, generation by 
convergence is more than offset by the effects of friction 
2 vc residual in Fig. 5.3 a) and eddy transport losses. 
The initial surface convergence in Fig. 5.3 d is onlya 
Mmeverron of that in Fig. 3.10 d for the developing case. 
This weak low-level inflow is consistent with the feeble 
vertical mass circulation (see Fig. 5.2 d-f). 

The strong NVA aloft at early stages of this non- 
developing case has no counterpart in the developing storms. 
Weak PVA is evident at early stages of both developing Type 
Mons (see Figs. 3.10 d and 3.13 d). The similarity of 
Ehe time series of quasi-Lagrangian vorticity tendency and 
Бист Су advection aloft (Fig. 5.4) suggests that vorticity 
changes are primarily the result of advection. aes ni 
cates a lack of mu tu sori icity generation. Thüs. 
vorticity changes aloft occur as the budget volume migrates 


to an area of different upper-level vorticity. Tl. - 
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suggests a lack of coupling between the surface disturbance 
(which is used to locate the budget volume) and upper-level 
features. Vertical coupling occurs via the secondary mass 
circulation when changes in the upper-level vorticity struc- 
Eure (increasing PVA aloft) occur in conjunction with 
changes (increases) in the strength of the vertical mass 
circulation. The juxtaposition of lower and upper features 
in this case is unfavorable for such coupling. The continued 
location of the surface disturbance beneath the upper-level 
convergence quadrant of the primary 300 mb jet is not condu- 
cive to vigorous mass circulation increases, such as 


occurred in the developing cases. 
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The FGGE storm also formed beneath the upper-level 
convergence (Fig. 3.14 d) associated with the equatorward 
exit region of the jet streak aloft (Fig. 3.6 a). Howeveus 
Strong vertical coupling in that case was established as the 
jet streak rapidly migrated to a position southwest of the 
surface disturbance (£16 7005 2 The vigorous mass circu- 
lation increases (Fig. 3.8) coincided vith the incurs4ou5225 
the divergent exit region of the jet into the budget volume 
(Cook, 1983). The presence of a downstream ridge (Fig. 3.4 
d) contributed to rapid growth Tor mE F In the ems 
developing case, the jet remains well to the north (Fig. 5.2 
a-c) and the vertical mass circulation remains weak (Fig. 
5.2 d-f). The maintenance of cyclonic vorticity aloft Lom 
east appears to inhibit the development of PVA. 

The heat budget shows large temperature increases 
throughout the troposphere (Fig. 5.5 a,b) that are almost 
entirely due to the horizontal eddy transport term (thermal 
acyeci ion) in Fates e oce Budget terms associated with 
the vertical mass circulation (diabatic heating and adas 
batic cooling) in Fig. 5.5 d-£f are small. Tl :s xs con ma 
with the lack of thermal vave amplification in Fig. 5.1995 
Thus, warming is primarily the result of the Southeast moves 
ment of the surface disturbance to a warmer location lc 
Riga. 5:13: ye By contrast, the leading terms in the heat 
budget for the developing Type I NOGAPS polar low (Fig. 
3.15) are diabatic heating and adiabatic cooling: This E 
cates that temperature changes are not entirely due to 
advective effects, which is consistent with the ampli ieii 
tion of thermal features in that case (517 рин 

The heat budget residual (Fig. 5.5 a-c) is remarks 
ably small at all stages of this non-developing d1Sturbanaen 
This attests to the accuracy of the numerical schemes Usem 
to compute the leading terms of horizontal eddy flux (vies 


includes the budget volume translation "veloces and 
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quasi-Lagrangian potential temperature tendency. The much 
larger residual for the developing case (Fig. 3.15 az€) 
arises from uncertainty introduced during the Vert m 
Interpolation woe temperatures from sigma to pressure 
surfaces. This residual is correlated with the net mean flux 
convergence (Fig. 3.15 d-f). Since this term is the Smee 
difference between the horizontal and vertical mean flux 


terms (which separately amount to several hundred degrees 


per day), this temperature error is magnified considerably. 
In comparison, the mean mass circulation for the non- 
developing case is small, vhich Tresults in smaller 


uncertainty. 

These budget analyses indicate that vorticity and 
thermal changes are largely the result of translation of the 
budget volume relative to features of the vorticity and 
thermal fields. The budget terms responsible for in-situ 
geopotential and thermal vave amplification are compara- 
tively small and are confined to lover levels. This Ra 
consistent with the weak development of the upper trough 
(Fig. 5.1 e,f). The analysis shows that this non-devel 
disturbance is shallow and only veakly coupled to the flov 
alo” 

The thickness tendency fields for this case (Fig. 


5.6) are quite different from those of the developing cases 


(Pues. oa In the developing cases, thermal 
advection contributes warming (cooling) ahead of (behind) 
the developing wave. As the low-level circulation becomes 
more intense, the amplitude of this dipole increases. 


Because adiabatic cooling and diabatic heating “ten 
balance to a high degree, net thickness changes in ки. 


developing cases occur which are approximately in phase with 


the advective contributions. Such changes are assoc NE 
with geopotential increases (decreases) in the upper ridge 
(Lrongh) which indicate self-amplification of the "up 
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wave. Thus, lower and upper-level features are strongly 


coupled in the developing cases. In contrast, the thermal 
advection patterns for the non-developing case (Fig. 5.6 
a-c) have no strong dipole structure. Rather, a general 


Lung trend at the initial time (Fig. 5.6 j) results from 
translation of the disturbance toward a warmer environment. 
B- oiabatic cooling (Fig. 5.6 d) and diabatic heating 
(Fig. 5.6 g) are small, in accordance with the heat budget 
157 5.50 d). This is consistent vith the veak secondary 
20 Circulation (Fig. 5.2 d-f). The resulting surface 
convergence does not produce the lov-level circulation 
Meem@eases required fora strong thermal advection dipole. 
Only at 204 h is there any hint of a self-amplification 
pattern (Fig. 5.6 k). However, the dipole amplitude is less 
than 100 gpm per day and the resulting amplification of the 
upper wave is slight (Fig. 5.1 e,f). 

This analysis suggests that a characteristic of the 
developing polar low which is absent for the non-developing 
case is the diagnosis of self-development, which is indi- 
cated by patterns of net warming (cooling) ahead of (behind) 
the low center throughout a deep layer. Self-development 
indicates a cooperative interaction between the growing 
Surface circulation and the amplifying upper-level features. 
In the non-developing case, thickness changes which indicate 
self-development are not present. This is consistent with 
the weak vertical mass circulation and indicates a lack of 
vertical coupling in the non-developing case. | 

Self-development is a manifestation of baroclinic 
development, 4000007 Отел ти regions of lov static 
5 lily, enhanced baroclinity and high relative humidity. 
Therefore, a key to understanding differences between devel- 
oping and non developing disturbances may lie in diagnosing 
these conditions. cons Tderation of these three criteria, 
an attempt was made to discriminate between developing and 


non-developing disturbances at early stages of formation. 
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At the incipient stage, Fie ere stabilıty index 
for the developing NOGAPS Type I case was &.7°K per 100 mb, 
which is close to the moist adiabatic value appropriate to 
the surface air temperature. For the non-developing distur- 
bance, the index increased from &.5 to 5.2°K per 100 mb 
aimee the first 12 h. The moist adiabatic rate for the 
same period increased from 4.6°K to 5.2°K per 100 mb as a 
result of the increasing surface temperature. Thus, there 
is little difference between the static stabilities of the 
developing and the non-developing cases. In both cases, the 
stability is close to the saturated adiabatic value. 

Another consideration is the initial baroclinity in 
the region of the two disturbances. A surface temperature 
difference of 14°C occurs between points 5? lat. to the 
north and 5? to the south of the non-developing disturbance 
center (Fig. 5.l a). For the developing case (Fig. 3.3 a), 
this difference is 12°C. Another measure of the baroclinity . 
is the vertical wind shear. The wind speed increases from 
£ han 5 m s-' at the surface to nearly 50 m s-* at 300 
mb for the non-developing case. The corresponding upper- 
level vind maximum for the developing polar lov is only 35 
ub (Fig. 3.2 a). Thus, 5077707:oclınıty of the non- 
developing disturbance is slightly stronger than for the 
developing case. 

Another important property of the initial state is 
the ambient relative humidity. The largest latent heat 
release in ascent regions of extratropical cyclones would be 
expected to occur when the low-level air is nearly satu- 
rated. At the initial times, the relative humidity in the 
75 "ty of the non-developing case (Fig. 5.7 d) is around 
70 5 onpared vith 80Z for the developing storm (Fig. 5.7 a). 
At later stages, the relative humidity of the developing 
25” approaches about 90% (Fig. 5.7 b) while that for the 


non-developing case decreases (Fig. 5.7 e). This difference 
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is because significant precipitation (not shown) occurs ne 
the developing storm whereas none occurs with the non- 
developing case. Thus, the near-saturation in Fig. 5./ b is 
the result of the condensation occurring in the ascending 
air: 

This analysis appears to suggest that the initial 
conditions for the non-developing case vere also favorable 
for baroclinic development, as in the developing case. That 
¡EN large baroclinity, high" hümudrr” and low stature 
Stability were present at the initial stages of the Type I 
non-developing polar low. This suggests that other factors 
are involved in the development of the Type I polar low. The 
location of the non-developing polar low beneath the conver- 
gence quadrant of the upper-level jet appears to inhibit the 
growth of the vertical mass circulation. This is consista 
with the lack of coupling between upper and lower features 


discussed earlier. 


B.. NON-=-DEVELOPING IPE IL DISTURB NCE 


ile Synoptic Description 


The storm described here is a slowly decaying 
surface vorticity maximum (Fig. 5.8 a-c) that is embedded 
within a strong northerly flow on the western periphery of a 
deep cyclone. This case has many of the characteristic gm 
the Type II lows described above, and is quite similar to 
observed polar lows that form as cyclonic vorticity Topec Sea 
the poleward side of the parent low. The vorticity cente EEN 
this case had earlier been the primary surface low center, 
but it had decayed aS a new, more vigorous low developed to 
the southeast. It is of interest that this new low R 
mature stage of the Type II storm described earlier (compare 
Figs. 4 lece Sando. 3 sane Polar lows occasionally develi 


from the old vorticity centers associated vrtU7:486- 77 
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occlodedesysteomsa Given the rather similar location a 
initial conditions as for the developing Type II case, the 
purpose of this section is to determine why this center did 
not develop while it was in the polar air stream. 

The non-developing vorticity center migrates cyclon- 
ically about the parent low and is eventually absorbed into 
the mean. flow (Fige usc At the start of this Gages 
surface vorticity values in excess of 8 x 1077370000 
present, with surface vinds up to 20 m/s near the center. 
The surface vorticity center is located beneath a broad 500 
mb trough (Fig. 5.8 d). A fast-moving short-vave trouguluud 
an associated 300 mb vind speed maximum (Fig. 5.9 a) are 
located about 700 km to the west of the surface center. 
After 24 h, the upper-level vorticity maximum has moved to a 
position southeast of the surface center (Fig. 5.8 e). This 
upper-level short-vave trough appears to be coupled with a 
rapidly developing surface trough that forms 1000 km Sou 
of the center (Fig. 15628) and translates cyclonicanmmm 
around the parent low ahead of the decaying center described 
here. This trough to the south is thus another example of a 
developing Type II polar low which forms in response to PVA 
aloft. Although the northern (dissipating) vorticity ce€n2 
lies poleward of the 300 mb jet (Fig. 5.9 a-c), a favorable 


coupling with a mid-tropospheric PVA center does not ocr okai 


Some upper-level support is initially present, as a center 
of upward vertical motion (Fig. 5.9 d) is located over the 
Surface center. However, this vertical motion weakens with 


time (Pier S orenk 

Very little baroclinity is associated with the mm 
ening disturbance (Fig. 5.10). In contrast, the Kora. EN 
temperature gradients associated with both the developing 
Type II disturbances (Figs. 4.3 and 3.3) are larger. espe 
cally at early =staces- The enhanced  baroclinity (07702 


SONS thu Tn ice 5.10 a is associated with the developing 
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disturbance there. The surface center at the initial time is 
located in the thermal ridge ahead of a thermal trough (Fig. 
5.10 a,d). Thus, the thermal patterns at the early stages of 
this non-developing disturbance are similar to those at the 
final stages of both Type I and Type II developments (Figs. 
EENM- t. 4.3956 .£:953.6 c,f and 4.5 c,f). This is consistent 
with the previous history of this non-developing disturbance 
as an occluded primary low center. 

Another indication G iO baroclinity 1s the 
vertical wind shear. The 300 mb wind speed directly above 
ee initial surface vorticity center is less than 10 ms-! 
(Fig. 5.9 a-c). The wind speed maximum in the jet does not 
exceed 40 ms-!. By contrast, wind speed maxima in the upper 
troposphere exceeded 55 ms-! above both the developing Type 
uubdisturbances (Figs. 4.2 a-c and 4.6 a-c). 

The surface circulation of this decaying storm is 
initially large. However, the lack of favorable coupling 
with the upper-level jet and weak baroclinity appear to 


inhibit further development of this disturbance. 


2. Quasi-Lagrangian Budget Analysis 


The quasi-Lagrangian budgets of mass, vorticity and 
potential temperature were calculated to determine the phys- 


ical processes that led to decay of this system. Although 


the horizontal analyses depict fields to 240 h, the 
quasi-Lagrangian analyses extend only to 216 h. After 216 
wend maximum of surface vorticity is no longer available to 


position the budget volume. 

The mass budget illustrates a layer of slowly 
decreasing convergence (Fig. 5.11 a) beneath a weakening 
outflow above about 700 mb. The deep layer of convergence 
Acne initial time reflects this disturbance s previous 
Bury as a major circulation center. The weakening 


2” I Velocity (Fig. 5.11 b) and the decreasing elevation 
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of the LND reflects the dissipating nature of this distur- 
bance and is consistent vith the significant but veakening 
upward motion in Fig. 5.9 d-f. Positive values of omega at 
the surface are consistent with increasing surface pressure 
that is associated with the vorticity center. Io contrast, 
the developing cases exhibit a rapidly growing mass circula- 
tion which sustains a deepening and strengthening convergent 
layer (Figs. 4.7 and 4.11) and spinup of vorticity. In the 
developing cases, the increasing ascent is coupled vith 
increasing PVA aloft near the poleward exit region of a 
Significant upper-level jet. Sc mm ame does not occur 
in the non-developing caSe Since the area of largest PVA 
aloft is coupled to the developing disturbance to the south. 


The vorticity budget is presented in Fig. 5.12. The 


SETONE, pre-existing vertical mass circulation tends to 
maintain a large, but slowly decreasing, convergence source 
term (Fig. D vci) near the surface. This low-level 


EN crty source is opposed by large frictional dissipation, 
which is implied by the large negative budget residual (Fig. 
7:7 :50). The friction is large because of the combination 
wenen surface winds and large positive vorticity. Even 
though the existing low-level convergence tends to sustain 
the vortex, the low-level vorticity decreases slowly (Fig. 
2— :)? In comparison, the convergence source for the 
developing Type II case increases rapidly (Fig. 4.9 d-f) and 
surface friction becomes large only at later stages of 
development (Fig. 4.9 d-f). 

ivemeasime voutieity aloft (Fig. 2) is the 
2” or PVA which is partially offset by divergence (Fig. 
E c.d). Docs “777 ni ii larger than that which 
accompanies the formation of the developing Type II storm 
Meer. 2 d). This large PVA would seem to be favorable for 
low-level cyclogenesis. However, an examination of Fig. 5.8 


2770 .-dicates that this large PVA is actually associated 
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with the migrating short-wave ии that is strong 
coupled with the developing disturbance to the south. In 
fact NVA occurs directly above the non-developing center. 
Thus, the seemingly favorable features of cyclogenesis and 
PVA aloft in the vorticity budget are associated with the 
developing system to the south that is also incorporated 
into the 5° latitude radius budget volume. 

The leading terms in the heat budget are diabatic 
heating and adiabatic cooling (Fig. 5.13 c,d). These teu. 
are about half those from the developing NOGAPS Type II 
Storm (Fie: 4.14 d-f), which reflects the larger vertical 
mass circulation for the developing case. Model cumulus 
heating (Fig. 5.13 e,f) is negligible for this case 102 
may indicate a model deficiency, Since in cold ponki 
outbreaks similar to Fig. 5.8 a,b, cumulus convect tome 
common. It may be that the type of convective overturn 
that occurs in nature also occurs here, but because MES 
confined to the lowest model layer it is labelled in NOGAPS 
as large-scale precipitation (e.g., Bosse, 1984). The NOTE 
zontal eddy term (thermal advection) is initially negative 
(Fig. 5.13 e). This is consistent with the cold advection to 
the south of the disturbance at low levels ,(Fig. 5.10 a) and 
the cold advection aloft that is associated with the short- 
wave to the west of the surface center (Fig. 5 M At 
later stages, the horizontal eddy term is smaller (Fie mE m= 
f) due to the more rapid southeastvard movement of the 
disturbance (see Fig. 5.8 a,b). Another difference betveen 
the heat budgets at early and later stages of this distur- 
bance is the decrease in the net mean flux convergence 
(compare Figs. 5.13 c,d). This decrease arises from 
increasing contribution from subsidence warming to tie Som 
of the disturbance center (Fig. 5.9 e) which is partia mE: 
included in ceune?budeee olum” This subsidence is associ- 


ated with the developing disturbance to the south. 1405770777 
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2 66 trom met cooling (Fig. 5.13 a) to warming (Fig. 5.13 
b) is due to increases in both the horizontal eddy and the 
net mean heat flux terms. As was the case with the vorticity 
budget, features associated with the developing storm are 
incorporated into the budget volume. This suggests that a 
smaller budget volume is required to better resolve the 
physical processes that are associated only with the dissi- 
pating vorticity center. 

The time section of diabatic heating (Fig. 5.14 a) 
shows an initial increase in heating followed by a dimin- 
ishing layer of heating which is strongest near the surface. 
Large surface fluxes that result from strong northerly flow 
over warmer water account for near-surface heating rates of 
around 4°C per day. For the developing Type II distur- 
bances, it was hypothesized that large surface heating 
combined with cooling aloft could enhance development of the 
upper trough. In this dissipating case, the upper-level 
550) is not ccupled vith the surface vorticity center, so 
such a coopera: ion between the surface heating and the 
amplifying upper-level trough cannot apply. 

meae n of baroclinic conversion (Fig. 5.14 b) 
is rather surprising in that the maximum conversion does not 
Eus until 204 h. A small positive conversion would be 
expected for the dissipating disturbance from the correla- 
wo ascent (Fig. 5.9 d) with the thermal ridge (Fig. 
5.10 a,c). The magnitude of both the thermal wave and the 
BO ucal circulation decrease at later stages seem to indi- 
cate reduced conversion. Thus, the increasing conversion in 
Ми 14 D during the first 12 h is probably due to the 
HüÜ0poration of the subsidence (Fig. 2775)? Wii eh occurs 
275: the amplifying cold trough (Fig. 5.10 b,e) to the 
South. This is coupled with the developing southerly 
Een. Thus, the large baroclinic conversion in Fig. 5.14 
b is due to the large budget volume that includes more than 


27 . “he circulation of the non-developing disturbance. 
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The horizontal distributions of terms in che emka 
ness tendency equation are shovn in Fig. 5.15, so thacozğ 
contributions from the tvo systems can be separated. The 
thermal advection dipole (Fig. 5.15 a-c) associated vıreudzur 
dissipating vorticity center is weak and is of {the wrens 
phase for self-development. The contribution from thermal 
advection is small because of the reduced baroclinity near 
the non-developing center. The net tendency (Fig. 5 > P 
results in increasing (decreasing) thicknesses behind (ahead 
of) the dissipating surface center. The development of a 500 
mb ridge to the northwest of the center (Fig. 5.8 e) k m 
is consistent with these changes, appears to augment the NVA 
directly above the center. These changes oppose the vertical 
mass circulation and development of the low-level trough. 
Ti contrast, thickness changes associated with the devel- 
oping Type Ll disturbance, (hie 4.18) are consistent with 
self-development and the growth of PVA aloft. 

The large thermal advection dipole (Fig. 5.15 22022 
and associated thickness tendency dipole (Fig. 5.15 5 0 5 m 
the south are linked to the developing disturbance there. 
The large dipole amplitude, which occurs because of the 
enhanced baroclinity in the south of the domain, is 247222 
able for self-development. As these patterns are partially 
incorporated into the 5° latitude radius budget volume gg 
vorticity and heat budgets centered on the dissipating 


system to the north must be carefully interpreted 
05 “Bisous-don 


The Type II polar low was shown to be a baroc mmm 
disturbance which is triggered by PVA aloft.  Developmen NEN 
associated with a developing upper-level trough Һати 
coupled to the surface cyclone. The non-developing Typemm 
disturbance studied in this section is a large-scale vores 


which decays as a result of Surfidee iriver 02 Although the 
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vorticity budget indicates large er 0 this is mop 
properly colocated with the dissipating cence: An ampli- 
fying short-wave trough aloft and associated PVA is coupled 
with another developing system to the south, which happens 
to be incorporated into the budget volume. The thickness 
changes associated with the decaying system actually oppose 
the development of upper-level PVA. Thus, the system slowly 
dissipates as the upper-level support decreases. As for the 
non-developing Type I disturbance, the absence of favorable 
coupling with the upper-level jet streak and associated PVA 


is the major factor which led to the decay of this system. 
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VI. SUMMARY AND CONCLUSIONS 


The purpose of this study has been to determine the 
reasons for the grovth of baroclinic disturbances vhich form 
vithin polar air streams. Formation and development of 
these polar lows have been studied using quasi-Lagrangian 
(storm-following) budgets of vorticity, heat and moisture. 
These diagnostic analyses have been applied to polar lows 
from two data sets. The first set is from the model simula- 
tions of Sandgathe (1981). These simulated polar lows formed 
Mene polar air stream of a mature extratropical cyclone. 
Therefore, the effects of nonlinear deformation processes 
associated with the parent lov are included at polar low 
formation. These conditions are more realistic than in the 
many modeling studies in which the disturbance forms on a 
zonal (unperturbed) basic state. The Second Set is the 
European Center for Medium Range Weather Forecasts (ECMWF) 
analyses of observational data from the First GARP Global 
Experiment (FGGE). This study extends the diagnostic anal- 
ysis of a polar lov previously studied by Cook (1983) and 
Wash and Cook (1985). A second polar low has been extracted 
from the FGGE analyses for in-depth studies. Since the FGGE 
analyses incorporate extensive supplemental observations 
over the oceans, they are especially suitable for the study 
of maritime veather systems. 

255: P- öy uses a combination of simulated and observa- 
wnal data. De ributiono: observational “data over 
the ocean is sometimes inadequate to resolve small-scale 
dısturbances. By contrast, model output forms a dynamically 
consistent data set vhich is available at every gridpoint to 
computer Troundoff accuracy at every time step. Such a 


Enpdete specification of physical processes cannot be 


Go 


achieved in an observational study. The increased time ии 
lution of model output was very beneficial in resolving the 
rapid changes in these polar disturbances. The completeness 
of the model-generated fields allowed evaluation of interpo- 
lation and other computational errors associated with the 
various diagnostic routines, and enabled closure of budgets 
to acceptably small residuals. While strong conclusions are 
possible only on the model budgets, the realism of the model 
simulations is demonstrated by simular although noisier, 
budgets from the observational data set. 

The use of cubic splines enables accurate estimates of 
instantaneous time derivatives without time averaging. 
Spatial averaging is avoided by the use of a bicubic spline 
fit to the horizontal data. Thus, these cubic spline £€22 
niques preserve the spatial and temporal resolution of the 
data, and reduce the truncation error associated with com 
lation of space and time derivatives. The improved space and 
time resolution associated with the use of cubic’ splines 
appear to be beneficial in the analysis of small-scale, 
Short-lived disturbances such as the polar low. 

Two types of polar low are identified from this studs ии 
the basis of the conditions existing during the formative 
Stage. The first (Type I) forms on the equatorward side of 


an upper-level jet streak under straight or anticyclonic 


flow aloft. The disturbance is initiated by low-level) DARAS 
clinic processes, and is characterized in its early Stace. 
by a weak, surface circulation with light winds. The Seco 


(Type II) forms in a strong surface polar air stresm0142 
cyclonic side of the upper-level wind maximum. This type 
initiated by positive vorticity advection E alofof 
Non-developing examples of both Type I and Type Il polar 
lows are also examined. Comparisons between developing and 
non-developing cases provide additional insight into reasons 


for polar lowantens 5155 


10: 


Бист ии arising from this study is that polar lows 
may form beneath straight flow aloft on the equatorward side 
of the jet. There is some controversy concerning the neces- 
5250010 a Short wave aloft during polar low. formation and 
intensification. The composite study of Mullen (1979) and 
Mmimyemcase Studies of Reed (1979), Locatelli et al., (1982), 
and Mullen (1983) showed that polar lows that form on the 
cyclonic side of the upper-level jet have a perturbed upper 
flow during early stages of development. These studies 
represent examples of Type II polar lows. However, the Type 
I polar lows identified by the present study form on the 
equatorvard side of the jet beneath straight flow aloft. 
Although a similar genesis structure has previously been 
Edo for extratropical cyclones  (e.g., Petterssen and 
Smebye's (1971) Type A cyclogenesis), and for so-called 
"medium-scale" disturbances vhich form east of Japan (Nitta 
27777”“mamoto, 1974; Chen et al., 1985), these conditions 
have not hitherto been applied to the polar lov. 

One of the main conclusions from this study is that both 
Type I and Type II polar lows form and grow primarily 
Pl basic baroclinic instability processes, and thus 
exhibit many of the features of larger maritime extratrop- 
ical cyclones. Genesis occurs in an area of enhanced baro- 
clinity that is associated with an upper-tropospheric jet 
streak. As the polar low intensifies, the baroclinity 
decreases. The vertical ascent for the model-generated cases 
exhibits a comma-shaped Signature. The comma tail is located 


in a surface pressure trough in an area of enhanced baro- 


IL. Another area of enhanced baroclinity occurs to the 
east of the surface lov. These features are similar to the 
£5  ."cal frontal cyclone. Similar frontal structure was 


Peewrously observed by Locatelli et al., (1982) for a number 


22 ih Pacific polar lovs. 
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This study shows that polar lows are associated M 
enhanced baroclinity and reduced static stability, whichis 
consistent with the conclusion that these disturbances are 
baroclinic phenomena. The small size of the polar low is 
due to the low Richardson number (indicated by the combina- 
tion of strong baroclinity and low staticuctacm mes as 
Suggested by Mullen (1979). Small-scale baroclinic dis CUMS 
bances are favored when the Richardson number is small near 
the surface (Gall, 1976a; Staley and 6211770 7772 

İT) SULOnHE” cold air outbreaks in which Type IL poki; 
lows grow, this study has shovn that enhanced surface fric- 
tion, and sensible heat fluxes may oppose formation of an 
unstable surface perturbation. Surface friction opposes the 
growth of low-level vorticity. Surface heating is largest in 
the coldest air, which tends to reduce the baroclinity. 
Addition of heat in antiphase with the growing temperature 
perturbation. associated with a baroclinically unscamm 
disturbance opposes development. Because of these stabi- 
lizing effects, it is suggested that the surface distür r m 
requires mass circulation increases which are triggered by 
PVA aloft. The stabilizing influence of the surface fluxes 
and the requirement for significant upper-level PVA possibly 
explains the rarity of polar lows compared to the frequency 
ol strong. Cold outbreaks. 

The composite study of Mullen (1979) indicates 
polar lows tend to form in areas of enhanced air-sea temper- 
ature difference, which implies a cyclogenetic role form 
surface fluxes. A hypothesis to explain this param 
regarding the role of the surface sensible heat flux in the 
case of the Type II polar lows has been suggested by this 
study. These polar lows tend to form in a low sta 
Stability environment. Strong surface heating may then lead 
to further destabilization throughout a deep layer. Ay on 


nation of strong surface heating and cooling alo tc 


cold sector could enhance development of the upper trough. 
Although strong surface fluxes oppose surface lov formation, 
it is hypothesized that strong surface heating in the cold 
air may indirectly contribute to polar low development via 
enhancement of the self-amplification process in the mid- 
troposphere. Type Ir polar lows form om the cyclonic side 
of the jet where a combination of reduced static stability 
through a deep layer, enhanced surface heating and PVA aloft 
create an environment conducive to development. con pa 
ison, the surface fluxes for the Type I disturbance are 
weak, and the initial perturbation, which is triggered by 
low-level baroclinic processes, will form near the maximum 
surface temperature gradient. Because an east-west baro- 
clinic zone slopes toward the pole, the Type I surface low 
will tend to form on the equatorward side of the upper jet. 
The budget analyses reveal several physical processes 
that are normally associated with a growing baroclinic wave. 
In both Type I and Type II polar lows, development is accom- 
panied by an increasing vertical mass circulation. The 
Mm@emeasine mass inflow at low levels results in vorticity 
27 ion, while increasing PVA aloft that is associated 
with an amplifying geopotential wave tends to offset the 
divergence sink. This upper-level wave is coupled with the 
intensifying surface low. The growth of the upper wave is 
more rapid for the Type II polar low. In the Type 1 case, 
generation of positive (negative) vorticity by convergence 
27  cCence) is less efficient in the lov absolute vorticity 
environment on the anticyclonic side of the upper jet. The 
budgets for the model-generated data are similar in most 
respects to those from the observational data. However, data 
uncertainty results in less reliable diagnoses of budget 
terms in the observational case. The vorticity budget 
25 716(Ss for both types are rather similar to studies of 


eer scale extratropical cyclones over or along the east 
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coasts of continents (e.g., Chen and Bosart, 1777779077” 
and Bosart, 1982, Calland, 1983, Bosart andi11 7 7 

The grovth of terms in the heat budgets of both types 
reflect the groving mass circulation. Diabatic processes, in 
particular condensation in the ascending air, augment the 
vertical mass circulation. This leads to increased vorticity 
generation near the surface and contributes to growth of the 
upper-level wave. À similar role for latent heat release has 
previously been described for extratropical cyclones (e.g., 
Smagorinsky, 1956; Aubert, 1957; Danard, 1965). Latent miea 
release in the warm sector of the developing polar low leads 
to an increase in the available potential energy (APE) of 
the system. Intensification of the polar low is accompanied 


by baroclinic conversion of this eddy APE to eddy kine ENS 


energy (KE), in accord with previous work for mid-latitude 
cyclones (e.g., Danard, 1966; Bullock and Johnson, ia 
Gall... 1976b:  Vincentrer als 1977: Robertson and Smir Tii 


1983). The numerical study of Hodur (1984) showed Piwi 
inclusion of condensation results in increased generation 
of eddy APE and conversion to eddy KE in Pacific polar ma 

A second conclusion of this study is that the sequens 
of self-development (Sutcliff and Forsdyke, 1950) tht 
been previously applied to larger scale exwtratropis mi 
cyclones (e.g., Palmen and Newton, 1969), applies to the 
developing polar lov. Self-development indicates a coopera- 
tive coupling betveen lover and upper-level features. 
Self-development is suggested by thickness change patterns 
of warming (cooling) ahead ot (behin the surface low 
center. Such changes are Mhhydrostatically consistent 7755 
amplification of the upper-level geopotential wave. 

In contrast with the developing cases, thickness change 
patterns in the non-developing cases studied here do not 
favor self-development. The lack of development in these 


cases is ascribed to the lack of coupling of lower NM 


upper-level features. The presence of a vigorous baroclinic 
disturbance just downstream of the non-developing system 
appears to oppose the mass circulation of the non-developing 
case by maintaining upper-level NVA and convergence. 

The budget results demonstrate the importance of the 
nonlinear effects due to thermal and vorticity advection. 
Mid-tropospheric vorticity advection is required for the 
genesis of the Type II polar low. Thus, the importance of 
nonlinear terms at polar low genesis cannot be dismissed. At 
later stages of development of both polar low types, strong 
deformation fields due to thermal and vorticity advection 
create conditions which appear to be favorable for future 
formation of secondary disturbances. Thermal advection 
contributes to frontogenesis in areas to the southwest and 
east, vhich are the expected locations of cold and varm 
F onts. These areas of enhanced baroclinity occur in 
conjunction with upper-level cyclogenesis. The associated 
PVA appears to initiate vertical ascent Which wou ld 
contribute to the development of a surface disturbance. It 
is suggested that diagnosis of frontogenesis and upper-level 
cyclogenesis may be the key to predicting the likelihood and 
50 I Ton of subsequent secondary low formation. 

This study indicates that the NOGAPS model does not 
appear to simulate the occlusion process vell. The observa- 
tional data exhibit a vell-defined peak in the strength of 
Wwe vertical mass circulation, with a decrease in the 
vertical slope of the geopotential features. By contrast 
the model appears to maintain the strength of the vertical 
Mass" circulation, and have only a modest decrease in the 
poucal siope. A second deficiency in the NOGAPS model is 
that the effects of cumulus heating appear to be underesti- 
used Cumulus processes result in the addition of heat 
higher in the troposphere than latent heat release arising 


217 large-scale ascent. Heating aloft and cooling near the 


171 


surface lead to stabilization of the atmosphere, which 
results in the demise of the storm. In the NOGAPS cases 
analysed in this study, large-scale ascent remains strong. 

Future studies of polar lows should focus on the condi- 
tions that lead to genesis. Because the polar low is a 
Short-lived phenomenon, a better understanding of genesis 
conditions would enable more timely forecasts of polar low 
formation. Further research should clarify the role 012557 
deformation fields of the parent low in the formation of 
secondary disturbances. Attention should also be focused on 
the evolution of low-level baroclinity in the presence of 
surface fluxes. Future modeling studies of polar lows should 
focus on the formation of secondary disturbances in pem 
airstreams to determine the role of processes which are 
associated with the parent low. Studies based on idealized 
zonal initial conditions may underestimate these nonlinear 
processes. 

Polar lows are poorly resolved by observational data. 
Data insertion schemes based on subjective interpretations 
of satellite imagery will be necessary in lieu of new data 
Sources. The vertical ascent signature provides information 
about the divergent component of the wind aloft. Objective 
analyses of polar lows should reproduce the vertical ascent 
Signature, and other satellite-viewed features such as cim 
direction of. cloud StEreeCs, cirrus striations and ci roe 
blow-off from cumulonimbus Fany a les Given these improved 
initial Conditions, the existing models appear to cona 
the physical processes for predicting the formation m 
intensification of polar lows of the type Studied “ewes 
Thus, there is some hope of improved forecasts of these 


Vigorous weather systems. 


j 


APPENDIX A 
DIAGNOSTIC FORMULATION 


A. QUASI-LAGRANGIAN BUDGET 


This study utilizes the geometry of the quasi-Lagrangian 
diagnostic budget framework formulated by Johnson and Downey 
1075). The budget equations are derived in a variation of 
spherical coordinates Lor nearly-cylindrical volumes 
consisting of Sections of the cone from the earth center 
contained between isobaric surfaces (Fig. А.1). The volume 
translates with velocity C, and moves with the storm. | 

The generalized isobaric budget expression of the quan- 
tity f is obtained by expanding the substantial derivative 
of f and integrating each term over the volume enclosed 
between pressure levels p and p, of horizontal area AA. 
This yields 
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Fig. A.l. Budget volume geometry. 


where f is the line-average of f around the circumference 
(AL) of the budget volume, £ is the area-average over the 
budget cross-sectional area (AA), and (V, -C,) is the invazo 
directed component of velocity normal to the trans 2727 
boundary. Budget terms are in units of f per unit true M 


deviation terms are computed as 


E cT E (at each circumference point) (A.2) 
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ет (at each grid-point) (A 


In (A.1), term (T) is the quasi-Lagrangian volume average 
Mate of change of f. Terms (HM) and (HE) are the gains of f 
due to convergence of mean and eddy horizontal fluxes. 
Terms (VM), (VE) represent increases in f due to vertical 
mean and eddy flux convergence (computed as the inflow 
through the lower boundary minus the outflow through the 
top). Term (S) represents the volume-averaged Sources and 
sinks, specified via the conservation equation for the quan- 
tity f. Term (R) is the budget residual. A positive residual 
ис that Ene observed increase of f is larger than is 


estimated by the sum of the computed terms. 


l. Computational Method 


goloulatıcno£ the cCerms in (A.T) follows Johnson 
and Dovney (1975). A generalized point P (see Fig. A.l) at 
polar coordinates (r,0) within the budget area has a lati- 
tude of 


.-l | l 
Ó = Sin İSin (6, 52) - (l-Sin 6)Sin (2) Cos D6! ии 


and longitude of 


Sin (2) Cos 9 


üs | (А.5) 


ANE AS + Sin 
O | Cos Y 


where @ and A are the latitude and longitude at the center 
B 7:57 7 the budget area, and a is the earth radius. All 
angles are in radians. Tnessnuinefe? Tr is Che arc length OP 
Methe earth surface, and the angle O is the azimuthal 
angle that this arc makes with the eastward-pointing vector 
77.) quanfity r/a in A.A and A.5 is the angle subtended 


at the earth center by the arc OP. Me component Jof 


B 


velocity normal to the boundary (positive ia dm at time 


ponmus OS s ivemibs 


V = ul 


T r,;9)Sin 8 - v(rg4,9)Cos B L 


where the angle f is 


E 7 Cosme 
б = Sin [Cos $6 Cos $. 


] (A.7) 
and the subscript B refers to a point on the boundary. The 
invard-directed component of C normal to the boundary is 


computed as 


"y 
(cəl (A.8) 
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с lel] сов(9 - tan " 


where C, and C, are the eastward and northward components oF 


C. The area average of f is calculated as 


r E 1 


f - [27a(1 ии Í i £(r,9)Sin(=)dédr (A.9) 
a 0 0 a 


The line average of f is 


ки Ада | 

ил. J Er ec (A.10) 
Simpson's 3/8 rule is used for both integrals. "Forms 
vertical integrals, a quadrature integration routine Dakki 


on a cubic spline representation of the integrand is used. 
This is more accurate than the simpler additive verre 
integration schemes commonly used. The circumference and 


area of the budget volume are computed as 


E 
as 2rasin (—) (A. 178 


and 


INS 


r 
AA = 2ma“ [1 - cos (2) ] ии 121 


To compute the terms in the budget expression (A.1) 
above, the meteorological variables are first interpolated 
from the rectangular lat/lon grid to the spherical storm- 
following budget grid. Equations (A.4) and (A.5) above are 
used to define the polar grid in terms of the rectangular 
lat/long grid. Then (A.1) is evaluated with the aid of equa- 
Moms (A.6) to (A.12). The budget grid used throughout is 
E umnensdoen 36 (amdimuth)- by 6 (radial) with a radial grid 
spacing of 1° latitude. 

Raw fields of NOGAPS output variables were found to 
be rather noisy, with apparently random grid-length scale 
perturbations superposed on the  synoptic-scale variations. 
Additional random uncertainty is introduced during the 
vertical interpolation from Sigma to pressure coordinates. 
zontal interpolation is performed by fitting a cubic to 


the data on the rectangular lat/long grid as described 


earlier. uan ts dara, this cubic will vary widely 
between data points, and intensify the noise on the budget 
erid. For these reasons, all model fields were smoothed 


EN to horizontal interpolation by using a 25-point 


smoother developed by Bleck (1965). The filter eliminates 
Short-wave noise while retaining the shape and magnitude of 
the synoptic-scale features. The filter also preserves the 
area average on the grid, an essential feature for a budget 
Emulation. 
2. Calculation of Kinematic Vertical Velocities 
mo ut us apparent that an estimate of 


75 (al velocity is required for the budget calculation. 
IhIS is obtained kinematically by vertical integration of 
the  divergences. Suitable boundary conditions must be 


imposed on omega at the top and bottom of the atmosphere. 
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For the NOGAPS model, omega is assumed zero at 50 mb. F am 
the FGGE analysis, omega is zero at 0 mb. Using the appro- 


priate boundary value, omega is computed at each point as 


P 
SC r ə | VAM dp (А.134) 


Pp 


A vertical spline fit to the vertical distribution of divas 
gence at each gridpoint is used to evaluate this integral. 
The technique avoids the cumulative error growth which 
occurs with the simpler additive downward integration 
schemes commonly used in the kinematic calculation. 

The bicubic spline fit to the u and v component Rm 
the wind is used to evaluate the divergences on the new 
grid. To ensure consistency with NOGAPS and FGĞGE formul. 
tions, earth curvature terms were added. The resulting 
divergences were adjusted slightly so that their area inte- 
eral was exactly equal to the line integral of the inflow 
velocity around the boundary. This eguali tye Eni Theorem), 
is a central assumption of the budget formulation. Prima 
adjustment, the two quantities differed by less th arms ии 
typical values. This excellent agreement attests to the 
mutual consistency of the interpolated wind and the diver- 
gence. Using the adjusted divergences, values of omega 
(vertical velocity) were then derived kinematically using 
(A.13) for each point om ehe”nev.erıü 

Because the model vinds are knovn to computer round- 
Ort waccupdgov the: “vertically integrated total colima 
divergences were used as the lower boundary value for omega. 
By comparison, it was considered that the substantial deriv- 
ative ot Surface pressure computed using the 
quasi-Lagrangian estimates of pressure tendency and transla- 
tion velocity from three-hourly data vas poorly p = 
Therefore, no adjustment was made to the divergences. Three 


estimates of the substantial derivative uvorelPlo:n9 MD 
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ib Fig. A.Z. The differences between these estimates 
reflect the effects of smoothing, interpolation, and time 
averaging. The instantaneous Eulerian pressure tendency is 
a model output field. The vertically integrated divergence 


yields realistic values of the substantial derivative, which 
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Fig. А.2. Three estimates of DP/Dt (mb h-!). 

Hozduue€ EGCGE data, the winds were considered to be 
poorly known. Therefore, the vertical velocities and diver- 
F ces were adjusted to fit both boundary conditions, using 
207 procedure outlined by O'Brien (1970). An analytic 


weighted correction to the vertical divergence profile was 


employed. This weighting function 


17” 


C(p) = (1000-P)/2000, “70 .8e. ии (A.14) 


has the form shown in Fig. Ann: 
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ô (Correction to divergence) 


Fig. A.3. Weighting function for correcting diversemeaam 


This function corrects preferentially near the surface) (game 
aos Above about /00 mb, the function increases linearly 
with height to allow for the increasing uncertaintey “527 
height of rawinsonde wind data. Wind observations near the 
surface contain much uncertainty as a result of differs 
anemometer exposures and boundary layer turbulence (wind 
gustiness) «far this reason: the strongest correction was 
applied below 900 mb, aS shown in Fie = SF The area- 
averaged divergences were corrected by an amount C(P) where 
Ö is obtained from 


= P 
ОР 5 


== + (V-€) vp, = f [9 u j ии A... 
öt S E 
Al, 
The second term in the integrand of AÀ.15 is Just ne rn 
gral of A.14. The corrected divergences were then used to 


compute omega using A.13. 
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The kinematic estimates of omega obtained by this 
method ensure mass balance at each point on the budget grid. 
This is necessary to avoid spurious sources and sinks in 
budgets of other quantities. In the vorticity, heat and 
moisture budgets, the horizontal and vertical transport 
terms associated with the mean circulation are each often an 
order of magnitude larger than the internal sources and 
sinks. However, their sum (net mean flux convergence) is 
rather small, and may be less than the sources and sinks. 
Therefore, in the present work, the mean mass circulation is 
balanced to better than 0.12 on the budget volume. iblis 
mass budget formulation differs from other budget studies in 
that mass balance is achieved at each gridpoint on the 
budget volume as vell as in an area-averaged sense. This 
Em ty is essential since mass continuity is an integral 
feature of model formulation at each gridpoint, and is 
assumed in the budget formulation. By this method, spurious 
Sinks and sources in the eddy terms arising from mass imba- 
lance are avoided. 

As a check on the reliability of the kinematic 
calculations for the NOGAPS fields, the kinematically 
derived omegas were compared with the model vertical motion 
ки 145. In sigma coordinates, 

dP 


Ç s 
= = + 5— . 
e (P, Prp) O Olay r xə 1 (5107 


7: € the first term on the right side of (A.16) is the 
model vertical motion in Sigma coordinates. The omega fields 
Calculated from (A.16) were interpolated directly to the 
pressure sufaces. Fields of omega (not shown) derived using 
(A.16) from the model vertical motion and computed using the 
kinematic method of this section are qualitatively very 
ar and differ quantitatively by only a fev percent. 


This close agreement suggests that the spline estimates of 
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horizontal derivatives are consistent with the finite 
difference approximation used in the NOGAPS formulation to 
compute the vertical velocity. However, to ensure mass 
balance, the kinematic rather than the model omegas are used 


for the budget calculations. 
jJ. Calculation of Translatron Ie e 


To position the budget volume, the maximum surface 
absolute vorticity is used initially as a first guess. The 
track obtained is then smoothed using.a least-Squares fit to 
a smooth cubic spline for both the time series of latitude 
and longitude. The smooth track so obtained departed by no 
more than half a degree latitude from the first guess. The 
Smoothed track positions are used to locate the budget 
volume. The x and y components of the instantaneous trans- 


lation velocity are obtained as 





C. = a cos ör, су a s (A.17) 
vhere and are the latitude and longitude of the budget 
volume at time t. The time derivatives are computed by 


differentiating the the cubic spline fit to the smoothed 
time-series of A and 0. This method produces a smoothly 


varying instantaneous translation velocity and mim na 


Spurious accelerations of the polar low. The instantaneous 
time derivatives of vorticity, heat, etc., are computed IR 
the same manner uSing a spline fit in time. These proce- 


dures enable computations of budget terms at each output 


time without time averaging. 
дад Transporte Partitioning 


The gains due to lateral and vertical transpor 
the budget quantity f are partitioned into mean and “0 


mode components, as shown in (A.l). The mean component is 


the change of f in the volume due to net flux convergence or 
divergence associated with the mean circulation. However, 
it is the small imbalance between the mean lateral inflow of 
259060 its vertical redistribution that contributes to the 
change of £. 

The horizontal eddy term represents the covariance 
of inflow and values of f on the perimenter. A positive 
eddy term results from inwardly directed wind components 
having large f values and/or outward flow with small f. The 


eddy term may be written as 





Vi ET < V f - Vn f (A.18) 























V = = Y. V (A.19) 
and 
| 
| m <= V- (EV) (A. 20) 
Therefore, 
VET = BP V.VE + SS(EV-V - f Vv] UNES 
Bus, the eddy term is equivalent to the area averaged 


advection of f only when the area average of the product of 
divergence and f equals the product of the area averaged 
divergence and the boundary mean value of f. This is approx- 
imately true when boundary and area averages of f agree 
closely and the perturbation quantities are small compared 
with the mean. Such conditions are best fulfilled with a 
210 budget volume. For a horizontal distribution of f 


and/or inflow that is symmetric with respect to the budget 


1: 


volume, the eddy term will be small. Thus, this term mi 
related to the asymmetry of the storm. The vertical eddy 
term represents changes in f due to the vertical difference 
between the spatial covariance of omega and f on the top and 
bottom surfaces of the budget volume. As such, it is quali- 


tatively similar to the vertical advection of f. 


B. CALCULATION OF QUASI-LAGRANGIAN ADVECTION TERMS 


Display of advection terms on a rectangular lat/long box 
poses a special problem in spherical coordinates because not 
all parts of such a box translate with the same velocity. 
For example, a box of fixed lat/long dimensions travelling 
north in the northern hemisphere shrinks in size. To remedy 
this problem, a correction is made to the transla 
velocity at each point within the box. (No such correct sen 
is required for the circular budget volume, since its arom 
is constant.) This correction can be substantial, espécr aie 
for a large display box having a large meridional trans" 
Trom velocity: The divergence in spherical coordinates may 


be written 





 — — - — tan 0 (A. 22 
Aa COSI. On 


Because every element of the box is translating with the 





constant vector velocity C - (с. c, ) it fol G 
dC RO 
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Geometrically, it is apparent that only the east-west dımana 
sion of the volume changes with north-south trans Cocose 
Therefore, the second term on the right hand side of MEE 


IS ZETO,UYL€C dr. 
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(A.24) 
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or 


C. (A) = C. + (À 5 Sin Фи (A025) 


where A-A, is the eastward displacement in radians from the 
meridian passing through the center point. 

Bü olculation of the vorticity advection term is 
described to illustrate the use of spherical coordinates. 


The relative vorticity in spherical coordinates is given by 
e E san 9 57 


il” x component of the vorticity gradient is obtained by 


differentiating (A.26) with respect to x to obtain 


2 
x Ly əl an 


ð dAd À 
E a 4 s 7 ф L 6660: 


O ÁÁ n 
— 


(A.27) 


Similarly, the y-component of the vorticity gradient is 





9C 2 ил E 
BB əz "95, 
7 a Cos o a dd a 
2555. (5225) 
GN OSEE 


7—— :.:):) using (4.25), (A.2/) and (A.28), the vorticity 


section is given by 





=. ад _ E 5—: 
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The derivative terms in (A.27) and (A.28) are obtained from 
the bicubic spline fit to the rectangular array of u and v 
data in ( , ) space and are mutually consistent. This tech- 
nique avoids the space-averaging inherent in finite differ- 
ence methods. The calculation of the thermal advection on 


the horizontal lat/long grid follovs7”” DF 
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